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Abstract 

The  functionality  of  modern  cities  relies  heavily  on  interdependent  infra¬ 
structure  systems  such  as  those  for  water,  power,  and  transportation.  Dis¬ 
ruptions  often  propagate  within  and  across  physical  infrastructure 
networks  and  result  in  catastrophic  consequences.  The  reaction  of  commu¬ 
nities  to  disasters  (e.g.,  seeking  alternatives)  may  further  transfer  and  ag¬ 
gravate  the  burden  on  surviving  infrastructures,  which  may  facilitate 
cascading  secondary  disruptions.  Hence,  a  holistic  analysis  framework 
that  integrates  infrastructure  interdependencies  and  community  behaviors 
is  needed  to  evaluate  a  city’s  vulnerability  to  disruptions  and  to  assess  the 
impact  of  a  disaster.  U.S.  Army  doctrine  requires  that  commanders  under¬ 
stand,  visualize,  and  describe  the  infrastructure  component  of  the  Joint 
Operating  Environment  to  accomplish  the  Army’s  missions  of  protecting, 
restoring,  and  developing  infrastructure.  To  this  end,  a  game-theoretical 
equilibrium  model  has  been  developed  in  a  multilayer  infrastructure  net¬ 
work,  to  systematically  investigate  the  mutual  influence  between  the  infra¬ 
structures  and  the  communities.  In  this  model,  two  types  of  infrastructure 
failure  patterns  are  formulated  to  capture  general  network  interdependen¬ 
cies;  network  equilibrium  is  extended  into  infrastructure  and  community 
systems  to  address  redistribution  of  demand  for  life-supporting  resources; 
the  societal  impact  of  disasters  is  estimated  based  on  resource  demand 
loss,  cost  increase,  and  total  infrastructure  failure.  A  real-world  case  study 
was  implemented  to  demonstrate  the  proposed  model  and  algorithm,  and 
to  reveal  insights. 
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1  Introduction 

1.1  Background 

Modern  cities  are  comprised  of  complex  infrastructure  networks  such  as 
those  for  power,  water,  transportation,  which  interact  with  one  another 
and  jointly  function  to  provide  resources  and  services  to  city  residents.  As 
cities  continue  to  expand  and  prosper,  the  ever-growing  population  im¬ 
poses  pressing  challenges  to  the  urban  infrastructure  systems  in  every  as¬ 
pect.  Even  for  properly  designed  infrastructures  that  satisfy  people’s  needs 
in  normal  functioning  scenarios,  infrastructure  performance  is  often  vul¬ 
nerable  to  unexpected  disruptions  due  to  factors  such  as  natural  disasters 
or  hostile  human  activities.  In  such  situations,  the  performance  of  the  city 
and  the  well-being  of  the  society  can  be  significantly  impacted,  resulting  in 
social  disruptions  related  to  economic  loss,  humanitarian  crisis,  and  de¬ 
mographic  loss. 

Many  events  in  recent  years  have  evidenced  that  infrastructure  failures, 
especially  those  under  cascading  propagations,  can  lead  to  catastrophic 
consequences.  A  well-known  example  is  the  Northeastern  blackout  in  the 
United  States  and  Canada  in  2003.  The  power  outage  started  in  Ohio  and 
then  propagated  to  several  nearby  states  and  the  province  of  Ontario,  Can¬ 
ada.  It  eventually  affected  over  50  million  people  and  led  to  a  total  esti¬ 
mated  economic  loss  of  $4.5  to  $8.2  billion  in  the  United  States  alone 
(Electricity  Consumers  Resource  Council  2004)  Thus,  the  impact  of  the 
blackout  spread  far  beyond  the  power  system.  First,  water  supply  systems 
were  impacted,  resulting  in  boil  water  advisories  across  many  regions 
(Johnson  and  Lefebvre  2003).  The  boil  advisories  were  mainly  due  to 
back-flows  in  the  water  treatment  systems  when  the  pumps  failed  to  pro¬ 
vide  sufficient  water  pressure  in  the  absence  of  a  power  supply.  Second,  in 
Ontario,  regional  gas  station  failures  occurred  due  to  malfunctioning 
pumps  and,  consequently,  many  long  queues  were  observed  at  the  remain¬ 
ing  working  gas  stations  and  fuel  prices  skyrocketed  (Public  Safety  2006). 
Third,  studies  also  showed  that  in  areas  such  as  New  York  City,  hospital 
admissions  and  healthcare  emergency  calls  increased  significantly  during 
the  blackout  due  to  many  kinds  of  induced  incidents  (Prezant  et  al.  2005 
and  Lin  et  al.  2011). 
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The  existence  of  complex  interdependencies  across  multiple  infrastructure 
systems  is  one  cause  of  the  fragility  in  urban  systems  under  disruption. 
Depending  on  the  geographical  location  and  function,  each  infrastructure’s 
functionality  relies  on  other  infrastructures  maintaining  their  functional¬ 
ity.  Particularly,  interdependencies  exist  among  components  (a)  within  the 
same  infrastructure  network  (e.g.,  water  treatment  plants  provide  re¬ 
sources  to  downstream  water  storage  facilities),  and  (b)  across  different 
infrastructure  networks  (e.g.,  power  plants  generate  electricity  for  water 
pumps,  traffic  lights,  and  fuel  refineries). 

Urban  infrastructure  failures  are  likely  to  stimulate  strong  reactions  from 
the  population.  A  direct  consequence  of  most  system  failure  is  difficulty  for 
residents  to  access  life-supporting  resources.  For  example,  people  may 
have  to  line  up  at  gas  stations  to  purchase  over-priced  fuel,  travel  a  longer 
distance  to  access  water,  or  turn  to  diesel  generators  when  the  power  grid 
is  disrupted. 

In  reality,  infrastructure  failure  and  community  reactions  are  mutually  de¬ 
pendent,  which  further  complicates  the  problem.  For  example,  people  may 
have  to  travel  through  the  transportation  network  to  deliver  or  retrieve  re¬ 
sources,  while  some  infrastructural  interdependencies  are  realized  by  de¬ 
livering  commodity  from  one  facility  to  another  via  transportation.  When 
congestion  increases  due  to  people’s  response  to  system  failure,  the  fluid¬ 
ity  of  commodity  flow  may  be  compromised  and  the  cascading  effect  could 
be  further  exacerbated.  Therefore,  instead  of  allowing  only  one-directional 
impacts  from  system  failure  to  population  response,  the  impacts  of  human 
activities  on  physical  system  performances  should  also  be  considered. 

Infrastructure  is  now  recognized  as  one  of  the  military’s  Joint  Operational 
Variables  and  Joint  doctrine-specified  missions  to  restore  essential  ser¬ 
vices  and  to  repair  and  protect  critical  infrastructure.  Furthermore,  in  the 
recently  released  “Army  Operating  Concept,  Win  in  a  Complex  World 
2020-2040,”  infrastructure  plays  a  key  role  in  the  descriptions  of  how 
landpower  is  leveraged  to  achieve  national  strategic  objectives  (U.S.  Army 
2014).  Doctrine  now  demands  that  commanders  and  staff  understand,  vis¬ 
ualize,  and  describe  the  infrastructure  operating  variable  to  accomplish 
missions  involving  selective  targeting,  protecting,  restoring,  and  develop¬ 
ing  infrastructure  in  combat,  stability,  humanitarian  assistance,  and  disas¬ 
ter  relief  operations.  Current  doctrine,  however,  does  not  say  how 
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commanders  and  staff  are  to  approach  these  challenging  tasks.  Current  in¬ 
frastructure  and  societal  models  do  not,  by  themselves,  provide  meaning¬ 
ful  representation  of  the  relationships  between  infrastructure  and 
populations.  This  shortfall  is  a  concern  for  both  military  and  civil  society, 
and  the  development  of  models  shall  work  equally  well  in  both  (Myers  et 
al.  in  process). 

1.2  Objective 

The  overall  objective  of  the  Human  Infrastructure  System  Assessment  for 
Military  Operations  (HISA)  research  package  is  to  develop  analysis  meth¬ 
ods  to  account  for  both  the  physical/functional  and  social  effects  of  infra¬ 
structure  changes  on  sub-populations. 

A  goal  is  to  build  a  model  that  represents  combined  human-infrastructure 
systems  so  that  the  potential  impacts  of  infrastructure  changes  on  social 
well-being  in  Army-relevant  contexts  can  be  explored.  This  model  will  be 
designed  to  provide  possible  policy  insights  into  how  best  to  protect  cru¬ 
cial  infrastructures,  reserve  emergency  supplies,  and  avoid  humanitarian 
disasters.  Example  applications  of  the  model  are  summarized  in  the  three 
following  scenarios: 

•  Destroy.  If  infrastructure  X  is  damaged  or  destroyed,  how  are  other  in¬ 
frastructures  affected  and  what  are  the  effects  on  end-users  in  the  com¬ 
munity?  Can  these  effects  be  represented  in  a  way  that  allows  negative 
impacts  on  local  populations  to  be  visualized  and,  as  a  result,  to  be 
managed  appropriately? 

•  Protect.  What  infrastructure  components  are  the  most  important  to 
preserve  or  defend  during  an  operation  if  the  well-being  of  a  particular 
community  is  to  be  prioritized? 

•  Restore.  Given  limitations  on  the  number  of  nodes  that  can  be  repaired 
or  rebuilt  in  a  given  time  frame,  how  can  repairs  or  rebuilds  be  priori¬ 
tized  to  have  the  highest  impact  (both  immediate  and  long-term)  on  lo¬ 
cal  communities? 

The  objective  of  the  subset  of  research  summarized  in  the  current  report  is 
to  develop  and  test  a  network  interdependency  model  that  provides  quan¬ 
titative  geospatial  representations  of  socioeconomic  impacts  of  changes  to 
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or  failures  within  an  infrastructure  system,  while  considering  that  popula¬ 
tion  reactions  to  infrastructure  failures  may  change  demand  patterns, 
which  themselves  affect  the  entire  system. 

1.3  Previous  work 

Complex  infrastructure  interdependencies  have  been  studied  by  a  number 
of  researchers.  Rinaldi  et  al.(200i)  discuss  the  challenge  of  modeling  and 
simulating  the  system  with  a  detailed  analysis  of  four  interdependency  cat¬ 
egories.  Zhang  and  Peeta  (2011)  review  work  addressing  infrastructural  in¬ 
terdependencies  with  different  approaches  and  propose  a  method  of  using 
computable  general  equilibrium  theory  to  generalize  the  interdependen¬ 
cies.  Focusing  on  different  types  of  interdependencies,  Motter  and  Lai 
(2002)  and  Buldyrev  et  al.  (2010)  respectively  performed  analytical  stud¬ 
ies  and  proved  that  disruptions  at  a  small  number  of  components  can  lead 
to  disastrous  cascading  failure  in  the  system.  Moreover,  efforts  have  been 
made  to  build  more  resilient  systems  (e.g.,  facility  location  and  design  cho¬ 
sen  to  optimize  reliable  service)  based  on  an  understanding  of  the  interde¬ 
pendencies.  For  example,  Wang  and  Ouyang  (2013)  describe  the 
geographical  interdependencies  of  demand  among  facilities  and  developed 
game-theoretic  models  to  find  an  optimal  spatial  distribution  of  service  fa¬ 
cilities;  Xie,  et  al.  (2015)  proposed  a  novel  network  augmentation  ap¬ 
proach  that  has  a  proven  ability  to  decompose  any  infrastructure 
disruption  interdependencies,  such  that  the  optimal  design  problem  can 
be  formulated  and  solved  equivalently  with  only  independent  disruptions. 

A  number  of  recent  studies  have  also  focused  on  the  social  impacts  caused 
by  people’s  behavior.  Horner  and  Widener  (2011)  first  considered  the  im¬ 
pact  of  a  transportation  network  failure  on  people  accessing  disaster  relief 
goods  after  a  hurricane,  and  then  tried  to  optimize  the  location  of  com¬ 
modity  distribution  facilities.  Many  efforts  have  been  made  to  model  com¬ 
munity  behaviors  following  an  emergency  or  a  disaster.  A  game-theoretical 
model  such  as  Nash  competition  has  been  used  to  capture  people’s  reac¬ 
tion  to  scarcity  of  resources.  In  particular,  it  is  probably  reasonable  to  as¬ 
sume  that  each  community  or  individual  would  try  to  acquire  resources 
with  the  least  time  and  monetary  cost,  taking  into  account  the  peer  reac¬ 
tions.  Such  a  situation,  often  referred  to  as  “equilibrium,”  changes  dramat¬ 
ically  after  infrastructure  failure,  and  those  changes  help  to  evaluate  the 
social  impact  of  that  failure.  Such  game-theoretical  models  have  been  ap¬ 
plied  in  related  literature,  e.g.,  within  the  scope  of  input-output  models 
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(Haimes  and  Jiang  2001;  Haimes  et  al.  2005).  Zhang  and  Peeta  (2011)  fur¬ 
ther  considered  an  equilibrium  that  took  into  account  both  spatial  and 
economic  factors  between  a  household  and  the  infrastructure  facilities  that 
provide  the  resources.  These  past  works  have  not  focused  specifically  on 
complex  interdependencies  of  multiple  infrastructure  systems,  and  this  re¬ 
port’s  authors  suggested  that  population  reaction  to  cascading  failures  be 
modeled  in  an  equilibrium  problem  by  considering  spatial  distributions  of 
facilities  and  population. 

1.4  Study  approach 

A  mathematical  framework  in  the  form  of  a  network  interdependency 
model  was  developed  to  systematically  quantify  the  effects  of  cascading 
failure  in  a  system  of  infrastructure  network  systems.  The  initial  approach 
was  developed  using  the  city  of  Maiduguri,  Nigeria,  as  a  case  study.  Using 
a  combination  of  existing  data  and  an  understanding  of  how  different  in¬ 
frastructure  components  are  used  in  this  location,  a  geospatially-refer- 
enced  network  model  of  the  critical  infrastructures  in  the  city  was 
developed  first.  Subsequently,  cascading  infrastructure  failures  were  prop¬ 
agated  under  different  scenarios  throughout  the  network  and  output  geo- 
spatially-referenced  before-and-after  costs  to  access  resources.  The  results 
of  different  failure  scenarios  were  then  assessed  and  several  types  of  sensi¬ 
tivity  analysis  were  performed  to  better  understand  both  system  and 
model  behavior. 

1.5  Study  location 

Maiduguri  is  the  capital  city  of  Borno  State  in  northeastern  Nigeria 
(ii°5i'N,  13°05'E),  with  an  estimated  total  population  of  1.2  million.  Con¬ 
current  with  rapid  urban  growth,  the  local  government  has  been  facing  ad¬ 
ditional  severe  challenges.  First,  the  natural  hazards  such  as  drought  and 
floods  cause  significant  adverse  effects  (Odihi  1996).  Next,  both  active  mil¬ 
itary  events  and  terrorist  attacks  threaten  people’s  daily  life  and  the  secu¬ 
rity  of  urban  infrastructure  (Ibeh  2015).  In  addition,  large  numbers  of 
internally  displaced  persons  (IDPs)  flee  into  Maiduguri  after  terrorist  at¬ 
tacks,  which  exhaust  the  resources  in  the  city,  resulting  in  further  pressure 
on  the  system  (Haruna  2015).  From  this  angle,  the  model  aims  to  better 
understand  and  interpret  these  pressing  social  concerns,  providing  possi¬ 
ble  policy  insights  into  how  best  to  protect  crucial  infrastructures,  reserve 
emergency  supplies,  and  avoid  humanitarian  disasters. 
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2  Methods 

2.1  Model  framework 

As  disruption  takes  place  in  an  infrastructure  system,  cascading  failure  be¬ 
gins  to  develop  by  following  certain  failure  propagation  patterns  that  are 
inherent  in  the  network,  and  people  start  to  change  the  way  they  access  re¬ 
sources.  Eventually,  the  disrupted  system  arrives  at  a  state  where  the  re¬ 
maining  infrastructures  and  the  impacted  population  behaviors  fall  into 
new  equilibria.  Modeling  and  querying  the  structure  of  these  new  equilib¬ 
ria  aims  to  understand  and  visualize  how  local  populations  might  be  af¬ 
fected  by  particular  infrastructure  changes. 

The  framework  therefore  (a)  generalizes  various  types  of  interdependen¬ 
cies  among  infrastructures  with  a  layered  network  model,  (b)  estimates 
entangled  system  failure  and  equilibrium  community  behavior  as  part  of  a 
game-theoretical  model,  and  (c)  evaluates  the  cascading  propagation  of 
disruptions  (due  to  interdependencies)  and  the  consequential  societal  im¬ 
pacts  (such  as  demand  loss  and  cost  increase).  A  heuristic  algorithm  is  de¬ 
veloped  to  calculate  the  system  equilibrium.  Numerical  comparisons  based 
on  real-world  data  are  performed  to  examine  the  impact  of  infrastructure 
cascading  failures  on  the  population,  and  model  sensitivities  are  explored. 

2.2  Network  development 

2.2.1  Nodes 

A  2-dimensional  geographic  region  O ,  such  as  a  small  country  or  a  city, 
was  the  focus.  Within  this  region,  the  infrastructure  system  is  character¬ 
ized  by  a  grand  network  that  contains  multiple  layers.  Each  of  these  layers 
represents  a  basic  type  of  infrastructure  (e.g.,  water,  food,  electricity,  fuel, 
transportation,  healthcare,  and  education).  It  was  assumed  that  each  layer 
consists  of  a  finite  number  of  infrastructure  nodes  as  the  functional  units, 
with  specific  geographic  locations  and  associated  attributes  for  those 
nodes.  For  example,  a  binary  “status”  attribute  at  each  node  describes 
whether  the  component  is  currently  functional  or  not.  Other  attributes  can 
be  included  as  necessary  to  further  identify  and  characterize  relevant  fea¬ 
tures  of  each  node,  including  the  nature  of  the  links  or  edges  that  connect 
the  node  to  other  components  within  the  system. 
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Figure  l  shows  four  example  infrastructure  sectors:  fuel,  power,  water,  and 
transportation.  These  sectors  serve  primarily  as  conceptual  boundaries 
that  help  analysts  define  the  scope  of  the  model  as  relevant  to  a  particular 
analysis.  The  sectors  are  not  fixed  or  hard-coded,  and  the  model  is  trans¬ 
ferable  to  locations  that  may  have  different  types  of  infrastructure  systems. 
In  addition,  each  sector  can  be  expanded  by  subdividing  it;  for  example 
the  transportation  sector  could  be  separated  further  into  air,  sea,  and  land 
systems.  Each  sector  is  composed  of  a  set  of  geographically-referenced 
nodes  representing  system  components.  In  this  example,  a  water  node 
might  be  a  water  tank,  a  fuel  node  might  be  a  gas  station,  and  a  transporta¬ 
tion  node  might  be  a  road  segment.  In  addition  to  the  four  layers  just  de¬ 
scribed,  a  community  layer  was  added  to  reflect  the  population 
distribution,  where  a  node  represents  a  functional  unit  appropriate  to  the 
scale  of  the  analysis  such  as  a  neighborhood,  household,  or  commercial 
entity.  Each  of  the  community  nodes  can  be  served  by  components  from 
the  infrastructure  layers. 


Figure  1.  Example  of  linked  infrastructure  nodes,  by  sector, 
with  community  end-users. 


Community 

Transportation 
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Sectors  such  as  food,  water,  fuel,  healthcare,  and  education  will  typically 
directly  provide  necessary  resources  to  the  community  layer,  hence  are  re¬ 
ferred  to  as  resource  layers.  Fuel,  electricity,  and  transportation  layers 
typically  serve  as  the  backbone  that  supports  the  operation  of  resource  lay¬ 
ers  in  addition  to  directly  supporting  community  nodes.  A  generic  node  in 
the  entire  network  is  denoted  by z  e  ^ .  Further,  let  r  e  R  be  a  resource  type 

(e.g.,  water,  electricity,  food,  healthcare)  and  ^  ^  be  the  corresponding 
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set  of  type-r  resource  nodes;  similarly,  let  ^  be  the  sets  of  transpor¬ 

tation  and  community  nodes,  respectively. 

2.2.2  Links 

After  the  major  functional  components  of  the  infrastructure  and  commu¬ 
nity  systems  have  been  defined,  links  can  then  be  established  between 
components,  either  within-sector  or  across  sectors  as  well  as  to  communi¬ 
ties.  In  order  to  account  for  both  direct  infrastructure  interdependencies 
as  well  as  system  feedbacks  induced  by  changes  in  community  behavior, 
the  model  currently  employs  two  different  types  of  linking  mechanisms 
representing  “strong”  and  “weak”  support  (Figure  2).  These  two  support 
types  then  determine  the  mechanisms  used  for  failure  propagation  in  the 
models,  described  in  later  sections. 

Figure  2.  Left:  In  strong  support,  if  one  of  the  three  required  inputs  to  the  top 
infrastructure  component  fails,  then  that  results  in  the  failure  of  the  downstream 
component  itself.  Right:  In  weak  support,  the  failure  of  one  resource  simply  means 
that  the  downstream  component  can  choose  from  among  alternative  sources,  and 
that  the  downstream  component  will  only  fail  if  the  access  costs  for  all  available 
sources  are  higher  than  some  predefined  threshold. 
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If  there  is  no  backup 

If  cost  is  too  much 

2.2.2. 1  “Strong”  support 

Where  the  functionality  of  a  “downstream”  node  has  a  direct  dependence 
on  an  “upstream”  component,  this  is  represented  in  the  network  model  as 
“strong  support,”  whereby  if  an  upstream  resource  fails,  the  failure  cas¬ 
cades  to  the  node  that  depends  on  that  resource.  There  currently  is  no  dis¬ 
tinction  between  functional  support  failure  (e.g.,  a  pump  cannot  move 
water  because  the  power  it  requires  to  run  has  failed)  and  input/output 
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disruption  (e.g.,  a  pump  cannot  move  water  because  there  is  no  water  to 
move  even  though  the  pump  itself  is  functional).  In  these  cases,  the  failure 
of  either  the  incoming  power  or  water  would  similarly  trigger  the  failure  of 
the  downstream  pump.  Typically,  strong  support  is  modeled  where  there 
are  direct  physical  linkages  between  infrastructure  components  that  neces¬ 
sarily  “push”  the  cascading  failures  downstream. 

2. 2. 2. 2  “Weak"  support 

Wherever  downstream  users  can  access  alternative  locations  of  a  resource 
via  the  transportation  network,  this  is  typically  represented  in  the  network 
model  as  “weak  support.”  For  example,  if  a  water  pump  has  failed  in  one 
location,  a  user  can  access  another  water  location  that  is  perhaps  farther 
away  but  still  relatively  convenient,  as  long  as  the  user  is  able  to  transport 
the  water  to  the  location  where  it  is  required.  In  weak  support,  down¬ 
stream  consumers  access  the  upstream  resource  based  on  a  user  equilib¬ 
rium  cost  calculation  over  the  existing  transportation  network.  However,  it 
is  noted  that  this  user  equilibrium  technique  is  not  restricted  to  transpor¬ 
tation  network  analysis,  and  can  represent  any  situation  where  resources 
are  accessed  via  some  network  capable  of  representing  those  costs.  Typi¬ 
cally,  weak  support  is  modeled  where  the  downstream  user  has  some  type 
of  choice  between  alternative  providers  of  a  resource  and  can  “pull”  the  re¬ 
source  from  the  most  cost-effective  site. 

2.3  Data  sources 

This  section  and  its  subsections  document  the  types  and  characteristics  of 
the  datasets  on  which  the  human-infrastructure  network  model  for  the 
Maiduguri  case  study  was  based.  The  mathematical  network  dependency 
model  with  the  solution  techniques  described  in  subsequent  sections  is 
general  and  theoretically  applicable  to  any  location  and  infrastructure  net¬ 
work.  However,  any  given  realization  of  the  model  will  rely  on  specific  data 
sources  to  represent  the  reality  on  the  ground,  and  the  utility  of  the  model 
will  depend  in  part  on  how  well  those  human-infrastructure  networks  are 
represented.  Given  that  the  modeling  tool  is  being  designed  to  be  applica¬ 
ble  to  areas  of  the  world  for  which  data  may  be  sparse  or  inaccessible,  in 
comparison  to  infrastructure  datasets  available  for  the  United  States,  it 
may  be  necessary  to  first  develop  the  input  datasets  from  imagery  and  in- 
situ  information  sources  in  an  area  of  interest  prior  to  model  formulation. 
In  future  work,  it  will  also  be  necessary  to  develop  methods  to  determine 
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the  conditions  under  which  model  users  can  have  confidence  that  the  in¬ 
put  data  available  to  them  is  adequate  for  the  analysis  of  a  given  situation. 

2.3.1  Urban  Tactical  Planner 

One  source  for  geospatial  data  on  infrastructure  in  urban  areas  outside  the 
United  States  is  the  Urban  Tactical  Planner  (UTP)  product  (U.S.  Army 
Corps  of  Engineers  2015).  This  dataset  provides  key  aspects  of  selected  ur¬ 
ban  environments,  including  built-up  terrain  zones,  buildings  of  interest, 
roads  and  other  infrastructure,  natural  features,  and  vertical  obstructions. 
The  primary  source  of  information  for  developing  the  UTP  product  is  aer¬ 
ial  imagery,  with  analysts  developing  several  vector  and  raster  geographic 
information  system  (GIS)  themes  from  the  image  base  and  other  sources. 
In  addition  to  roads  and  hydrography,  a  few  vector  files  typically  included 
in  UTP  which  may  be  unique  to  the  UTP  product  and  particularly  useful  to 
HISA  analysis  include  Key  Features,  Buildings  of  Interest,  Vertical  Ob¬ 
structions,  and  Built-Up  Terrain  Zones  (BTZs).  Each  of  these  layers  is  de¬ 
scribed  in  more  detail  below.  In  addition  to  these  vector  layers,  three 
raster  layers  are  included  in  every  UTP  product:  a  Base  Image,  a  Base 
Map,  and  an  Elevation  Layer.  For  some  locations,  the  UTP  product  may 
also  contain  layers  that  identify  where  different  ethnic,  religious,  or  politi¬ 
cal  groups  reside  (if  such  information  is  available). 

Since  UTP  products  are  developed  by  the  U.S.  Army  Geospatial  Center  to 
support  military  operations  in  urban  terrain,  UTP  products  only  exist  for 
select  urban  areas  where  existing  or  potential  conflict,  peace-keeping,  or 
humanitarian  missions  are  active  or  anticipated  from  the  perspective  of 
the  United  States  Department  of  Defense  (DoD).  Therefore,  an  UTP  prod¬ 
uct  does  not  exist  for  all  urban  areas.  The  UTP  data,  as  per  requirements  of 
this  project,  is  UNCLASS/For  Official  Use  Only  (FOUO),  with  access  re¬ 
stricted  to  DOD  common  access  card  (CAC)  users. 

Key  Features  is  a  point  feature  class  identifying  the  location  of  cultural 
features  (e.g.,  houses  of  worship)  or  operationally  significant  features  (e.g., 
police/government  buildings).  Of  particular  note,  the  Key  Features  layer  is 
one  of  the  major  sources  of  information  in  the  UTP,  describing  buildings 
and  facilities  that  are  protected  under  the  laws  of  war  (U.S.  Army  Corps  of 
Engineers  2013). 

Similarly,  Buildings  of  Interest  is  a  polygon  layer  that  shows  the  ground 
footprint  of  buildings  deemed  by  an  analyst  to  be  of  particular  importance. 
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Similar  to  the  Key  Features  point  layer  (and  often  representing  the  same 
structures),  the  Buildings  of  Interest  layer  delineates  high-profile  cultural 
features,  government  and  religious  buildings,  and  all  types  of  structures 
protected  under  the  laws  of  war  (U.S.  Army  Corps  of  Engineers  2013). 

Vertical  Obstruction  information  consists  of  vector  feature  classes  (points 
and  lines),  originally  of  primary  interest  to  pilots  for  the  purpose  of  avoid¬ 
ing  structures  that  could  impede  safe  takeoffs  and  landings.  These  features 
typically  include  infrastructure  components  such  as  radio  towers,  water 
towers,  and  power  lines;  therefore  the  layer  may  be  of  interest  for  HISA 
analysis.  These  layers  aim  to  include  all  point  structures  >15  m  in  height. 

Built-up  Terrain  Zones  are  delineated  by  using  a  method  of  grouping  ur¬ 
ban  terrain  according  to  its  use  and  physical  features.  BTZs  are  classified 
by  form,  function,  and  height. 

•  Form  is  categorized  as  Open-Spaced,  Detached  Closed  Space,  and  At¬ 
tached. 

•  Function  is  classified  as  Residential,  Commercial,  Manufacturing,  and 
Industrial  or  Institutional. 

•  Height  is  classified  into  height  ranges  of  1-2  stories,  3-4  stories,  5-15 
stories,  and  >15  stories. 

For  analysis  here,  the  function  of  BTZs  is  of  particular  interest.  The  resi¬ 
dential  function  type  includes  single-family  homes  and  apartment  com¬ 
plexes  along  with  temporary  housing  such  as  campgrounds  and  refugee 
camps.  In  some  cases,  the  residential  areas  are  labeled  according  to  socio¬ 
economic  status,  such  as  middle  income,  poor  income,  or  shantytowns. 
This  income  classification  is  based  upon  a  generalized  standard  of  income 
on  a  country-by  country  basis.  The  Commercial  function  type  identifies  ar¬ 
eas  that  are  primarily  used  for  economic  activities,  such  as  shopping  areas, 
markets,  office  complexes,  and  recreational  facilities.  It  also  includes  large 
transportation  complexes  and  areas  of  mixed  use  (e.g.  first-floor  commer¬ 
cial  with  apartments  above).  Manufacturing  and  Industrial  function  types 
include  large  factories  and  warehouses.  Institutional  areas  are  used  pri¬ 
marily  for  government,  religious  or  military  functions. 

2.3.2  Open  Street  Map 

One  open  source  for  global,  general,  background  geospatial  data,  including 
streets  and  oftentimes  hydrography,  electrical  lines,  and  railroads,  is 
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OpenStreetMap  (OSM)  (www.openstreetmap.org).  OSM  is  an  initiative  to 
create  and  provide  free  global  geographic  data.  Because  OSM  is  commu¬ 
nity  driven  or  crowd-sourced,  it  is  considered  a  non-authoritative  source 
for  geospatial  data.  However,  many  of  its  community  members  and  con¬ 
tributors  are  GIS  professionals,  and  in  some  situations  the  geospatial  data 
content  may  represent  best-available  data. 

2.4  Conceptual  model  construction 

This  section  describes  a  conceptual  model  of  the  infrastructure-commu¬ 
nity  system  in  Maiduguri  in  which  the  team’s  goal  was  to  find  a  balance 
among  an  accurate  depiction  of  reality,  best  use  of  available  data,  and  in¬ 
ference  of  additional  details  that  could,  in  an  operational  situation,  be 
more  robustly  supported  by  additional  in  situ  information.  Figure  3  shows 
the  conceptual  model  of  the  infrastructure  network  that  was  used  in  the 
subsequent  network  equilibrium  analysis  in  Maiduguri.  The  subsections 
below  describe  the  reasoning  behind  the  setup  of  the  various  sectors,  as 
well  as  highlighting  some  of  the  limitations  that  model  users  will  need  to 
consider  when  working  with  analyses  derived  from  this  model.  In  the 
HISA  workflow,  the  conceptual  model  serves  to: 

•  organize  analyst  understanding  of  the  infrastructure  system, 

•  guide  the  process  of  linking  geospatial  data  into  a  properly  constructed 
network  amenable  to  analysis, 

•  provide  a  reference  to  aid  in  the  interpretation  of  model  results,  and 

•  provide  a  record  of  how  the  analysts’  understandings  of  the  network 
evolved  over  time  as  revisions  were  made  and  new  information  was  in¬ 
corporated. 

The  conceptual  model  also  suggests  what  could,  in  the  future,  be  devel¬ 
oped  into  a  graphical  user  interface  (GUI)  whereby  the  creation  of  the 
mathematical  inputs  to  the  network  model  could  proceed  directly  via  user 
interaction  with  graphical  components. 

While  based  in  reality  as  per  the  data  sources  described  above  and  external 
research,  to  a  large  extent  the  conceptual  model  elaborated  for  Maiduguri 
is  a  simulated  dataset  that  shares  many  of  the  general  characteristics  of 
known  infrastructure  in  this  location.  However,  where  “proxy”  infrastruc¬ 
ture  components  are  identified,  this  means  that  the  actual  locations  are 
not  known  and  simplifications  have  been  made  to  arrive  at  this  particular 
component  configuration.  Each  sector  is  documented  below,  first  with  a 
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general  background  of  the  infrastructure  sector  followed  by  a  summary  of 
how  the  available  information  was  used  to  create  the  conceptual  model. 

In  the  conceptual  model  created  of  Maiduguri’s  infrastructure,  six  utility 
sectors  are  represented  (Fresh  Water,  Fuel,  Electricity,  Transportation, 
Food,  and  Community),  as  well  as  two  institutional  sectors  (Schools  and 
Medical  Facilities)  that  are  treated  as  infrastructure  components  rather 
than  community  end-users  (Figure  3).  The  sectors  in  the  figure  are  color- 
coded,  and  component  classes  within  each  sector  are  represented  by  col¬ 
ored  boxes.  These  boxes  are  then  interconnected  either  via  strong  support 
(e.g.,  “physical  links,”  with  solid  black  lines)  or  weak  support  (e.g.,  “trans¬ 
portation  links,”  with  dashed  black  lines)  to  other  component  classes.  Two 
end-user  communities  and  the  links  connecting  them  to  the  infrastructure 
network  are  also  represented  in  this  diagram;  other  end  users  were  repre¬ 
sented  when  the  network  was  encoded,  but  they  are  not  represented  here. 
Components  that  are  semitransparent  also  have  not  yet  been  included  in 
the  model.  The  additional  text  below  each  component  box  details  some  of 
the  attributes  that  characterize  the  linkage  types  for  each  component.  The 
following  section  describes  the  location-specific  model  representations 
and  assumptions  for  each  infrastructure  sector  represented  in  this  dia¬ 
gram. 


ERDC  TR-16-11 


14 


Figure  3.  The  conceptual  infrastructure  network  model  for  Maiduguri.  Colored  boxes 
represent  types  of  nodes  in  the  network;  solid  black  lines  and  arrows  indicate  strong 
support;  dashed  black  lines  indicate  weak  support.  Semi-transparent  nodes  have  not 
yet  been  represented  in  the  model  (ERDC-CERL). 
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2.4.1  Fuel  sector 

2.4.1.1  Background 

There  is  one  location  in  Maiduguri,  termed  the  “depot  terminal,”  where 
diesel  fuel  is  delivered  via  pipeline  to  the  outskirts  of  the  city.  The  source 
of  this  diesel  is  refineries  in  the  southwestern  part  of  the  country.  A  net¬ 
work  of  both  formal  and  informal  fuel  distribution  stations  exists  within 
the  city  to  distribute  this  resource,  many  of  which  may  be  co-located  with 
service  stations  identified  in  the  UTP.  Recent  news  reports  indicate  that 
Boko  Haram  insurgents  operating  in  Maiduguri  have  created  their  own  lo¬ 
cations  at  which  fuel  is  stockpiled  for  their  own  use  (Sahara  Reporters 
2015).  This  demonstrates  the  dynamic  nature  of  infrastructure  systems  in 
conflict  zones  and  underscores  the  importance  of  understanding  and  mon¬ 
itoring  both  formal  and  informal  infrastructure  networks  in  urban  areas 
where  hostile  and  friendly  forces  exist  in  close  proximity. 

2.4. 1.2  Conceptual  model 

The  conceptual  model  begins  with  the  inclusion  of  the  depot  terminal. 
Components  of  the  fuel  sector  external  to  the  city  itself  are  not  modeled, 
and  it  is  assumed  that  if  the  depot  terminal  is  operational,  then  diesel  fuel 
is  accessible.  A  subset  of  eight  locations  of  otherwise  unspecified  types  of 
“tanks”  in  the  UTP  that  were  located  near  roads  were  selected  as 
representative  locations  for  fuel  storage  and  distribution  in  the  city.  This 
first  iteration  has  not  integrated  service  stations  into  this  network  of 
possible  distribution  stations. 

The  UTP  does  not  include  information  on  the  layout  of  any  underground 
fuel  pipelines  within  the  city,  and  the  link  between  the  depot  terminal  and 
the  diesel  tanks  has  been  modeled  as  a  transportation  link,  presuming  that 
this  transfer  typically  takes  place  via  tanker  trucks.  It  is  important  to  note 
that,  currently,  the  model  does  not  track  actual  volumetric  flow  or  storage 
of  the  fuel  itself,  nor  volume  for  water  or  any  other  resource.  Because  each 
diesel  tank  relies  on  the  depot  terminal,  a  disruption  of  the  single  “depot 
terminal”  means  there  is  no  secondary  location  from  which  the  tanks 
might  pull  fuel.  In  this  case  therefore,  an  immediate  disruption  of  the  fuel 
supply  at  the  diesel  tanks  would  occur,  irrespective  of  how  much  fuel 
might  be  stored  in  the  tank  or  the  rate  at  which  it  may  be  depleted.  Diesel 
is  itself  a  resource  used  by  all  the  vehicles  in  the  transportation  system. 
However,  this  dependence  is  not  explicitly  modeled  at  this  time;  rather,  all 
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fuel  is  passed  via  weak  linkages  from  the  8  fuel  tanks  to  the  set  of  137 
proxy  diesel  generators  placed  in  commercial  locations  throughout  the 
city. 

2.4.2  Electricity  sector 

2. 4.2.1  Background 

Electrical  generation  occurs  in  the  southern  part  of  the  country  and  is 
transmitted  long  distances  via  330  or  132  kV  transmission  lines,  periodi¬ 
cally  boosted  at  transmission  stations.  Delivery  currently  comes  from  one 
transmission  line  that  goes  through  the  town  of  Damboa  to  the  south,  alt¬ 
hough  there  are  plans  for  an  additional  transmission  line  to  connect  from 
the  west  to  make  the  system  more  resilient.  Electricity  provision  typically 
does  not  keep  up  with  demand,  and  outages,  both  forced  and  planned,  are 
common  nationwide  throughout  the  network.  Electrical  infrastructure  has 
also  been  a  target  of  vandalism  and  concerted  attacks.  In  June  2014,  the 
330/132  kV  substation  in  Damboa  was  set  on  fire,  and  in  mid- July  2014,  a 
conductor  on  a  transmission  line  was  cut  somewhere  between  Damboa 
and  Maiduguri,  resulting  in  a  complete  electrical  blackout  for  Maiduguri. 
The  blackout  lasted  almost  three  months;  finally,  with  the  help  of  federal 
forces  for  security  and  logistics,  the  transmission  company  was  finally  able 
to  restore  power.  Distribution  is  managed  at  the  local  level,  wherein  Mai¬ 
duguri  is  one  of  four  districts  belonging  to  the  Yola  Electricity  Distribution 
Company.  After  various  privatization  initiatives  across  Nigeria,  Yola  now 
has  a  technical  partner  in  the  Philippines-based  MERALCO.  Yola  owns  a 
total  of  124  circuits  and  related  substations  (33  and  11  kV)  and  close  to 
1,000  smaller  distribution  substations.  The  UTP  identifies  a  set  of  pylons 
that  bring  electricity  from  the  high-voltage  country-level  transmission  net¬ 
work  to  the  two  nearly  co-located  electrical  substations  in  the  center  of 
Maiduguri  as  well  as  a  dual  substation  also  located  in  the  center  of  the  city. 
Sources  do  exist  with  names  and/or  descriptions  of  electricity  service 
zones  (e.g.,  Waziri  2009),  but  the  information  is  not  currently  available  in 
formats  amenable  to  geospatial  analysis. 

Several  internet-based  news  articles  provided  some  indication  of  the  ef¬ 
fects  of  the  summer  blackout  of  2014  as  well  as  other  infrastructure  at¬ 
tacks,  many  of  which  have  been  incurred  by  the  insurgent  group  Boko 
Haram. 
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Boko  Haram  has,  in  the  past  three  years,  vandalized  public  infrastructure 
like  telecommunication  masts  that  had  hitherto  cut  off  various  parts  of 
Borno  State  from  the  rest  of  the  world.  The  group  had  also  used  bombs 
and  fire  to  destroy  schools,  hospitals,  police  offices,  barracks  and  even 
cratering  of  roads.  Their  most  recent  attack  on  public  infrastructure  was 
the  bombing  of  an  ultra-modern  drilling  rig  procured  by  Borno  State  at 
the  cost  of  over  N300  million  at  the  site  where  it  was  mobilized  to  drill 
water  for  rural  dwellers  of  the  northern  part  of  Borno  State.  (Boko  Ha¬ 
ram  Attacks  Power  Station  2014) 

Simply  due  to  the  extended  length  of  the  Maiduguri  blackout,  it  is  likely 
that  private  operators  of  water  pumps  now  have  gained  gas-powered  gen¬ 
erator  access,  and  for  this  reason,  interdependency  between  the  electrical 
and  water  systems  may  be  quite  weak  compared  to  that  of  “more-devel¬ 
oped”  locations.  However,  another  researcher  notes  that  running  small 
businesses  on  generators  is  typically  prohibitively  expensive  and  that  the 
lack  of  electricity  at  night  “creates  a  cover  for  insurgents  to  operate”  (Mur¬ 
dock  2014).  One  article  also  highlighted  the  effects  of  the  blackouts  on 
subpopulations  that  the  model  does  not  explicitly  delineate:  “Cold  chain 
stores  for  fruits  and  vegetables  preservation  will  come  to  a  halt,  with  in¬ 
creasing  losses  of  goods  and  sales  revenues  of  marketers  and  traders,  in¬ 
cluding  the  over  15,000  butchers  and  Maiduguri  Monday  Market  traders 
that  solely  rely  on  trading  to  survive  and  send  their  children  to  school,” 
and  that  “the  economic  situation  of  vulcanisers  and  welders  that  consti¬ 
tutes  about  45  per  cent  of  labor  force  will  be  made  worse”  (Olugbode 
2014). 

2A.2.2  Conceptual  model 

The  conceptual  model  includes:  (a)  one  of  the  pylons  representing  the 
electricity  provided  by  the  national  network,  (b)  one  substation  to  repre¬ 
sent  the  two  nearly  co-located  nodes  in  the  UTP,  and  (c)  a  set  of  166  “proxy 
transformers”  throughout  the  city  in  institutional  and  residential  commu¬ 
nity  locations  to  represent  likely  community  links  to  the  power  grid.  The 
proxy  transformers  are  also  connected  to  schools  and  medical  facilities  via 
“strong  support”  links.  In  the  conceptual  model,  the  water  treatment 
plants  and  a  small  number  of  industrial  locations  are  also  provided  power 
via  the  substation  directly. 
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In  addition  to  the  fixed  electrical  infrastructure  network  that  is  notoriously 
unreliable,  many  commercial  enterprises  employ  diesel-powered  genera¬ 
tors  for  their  businesses.  The  conceptual  model  represents  this  by  placing 
a  generator  at  each  of  137  commercial  locations  identified  by  the  UTP 
across  the  city,  each  of  which  can  “pull”  diesel  fuel  from  one  of  the  8  diesel 
tanks  via  a  weak  support  link. 

2.4.3  Water  sector 

2.4.3. 1  Background 

Before  1985,  urban  access  to  water  in  Maiduguri  was  via  65  deep  boreholes 
with  pumps;  however,  two  severe  droughts  in  the  1970s  brought  rural  resi¬ 
dents  to  makeshift  urban  settlements  and  urban  groundwater  levels  de¬ 
clined  up  to  7  m  between  1975  and  1980  (World  Bank  1996).  This 
prompted  the  World  Bank  to  fund  a  project  costing  $72  million  in  U.S. 
dollars  (USD)  from  1985  to  1995  for  water  pumping  and  treatment  from 
Lake  Alau,  Nigeria,  located  southeast  of  Maiduguri.  A  water-sharing  agree¬ 
ment  specified  that  during  normal  years,  the  water  was  to  be  shared 
equally  between  urban  and  irrigation  uses,  and  during  drought  years,  pri¬ 
ority  would  be  given  to  urban  use.  The  project  included  a  treatment  plant, 
pumping  and  storage  facilities,  and  transmission  and  distribution  lines. 

Despite  completion  of  a  substantial  portion  of  the  planned  infrastructure 
improvements,  at  project  closeout  the  World  Bank  judged  the  outcome  to 
be  unsuccessful  for  a  variety  of  reasons  (World  Bank  1996).  These  reasons 
included  use  of  sensitive  equipment  ill-adapted  for  windy  and  hot  atmos¬ 
pheric  conditions,  and  possibly  seasonally-variable  raw  water  quality  con¬ 
ditions;  failure  to  establish  a  reliable  supply  chain  for  obtaining 
replacement  parts  and  chemicals  required  for  operation;  failure  to  procure 
resources  to  purchase  fuel,  spare  parts,  and  chemicals;  lack  of  monitoring 
and  control  of  the  water  distribution  network  with  failure  to  address  illegal 
connections;  unsustainable  costs  for  customers  (e.g.  sustainable  rates  to 
cover  utility  costs  would  have  required  charging  30%  of  an  average  civil 
servant’s  salary);  failure  to  train  management  and/or  consider  difficulties 
of  working  with  the  underperforming  civil  service  sector;  and  “too  much 
reliance  on  external  technical  assistance  to  perform  key  functions  of  pro¬ 
ject  implementation”  (World  Bank  1996)  Even  during  the  project,  one- 
third  of  boreholes  were  inoperative  due  to  lack  of  maintenance  and  only 
one  of  four  valves  at  treatment  plant  was  operable  at  time  of  project  close¬ 
out.  The  World  Bank  closeout  report  states  that 
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. .  .the  key  lesson  doubtless  remains  that  it  does  not  suffice  to  construct 
physical  facilities  for  provision  of  water;  of  equal  importance  -  if  not 
more  so  -  is  to  develop  the  institutional  capacity  of  the  water  utility  and 
to  provide  it  with  the  necessary  independence  to  successfully  operate  and 
maintain  physical  installations.  For  if  such  capacity  does  not  exist,  as  is 
the  case  with  the  Bomo  State  Water  Corporation,  expensive  installations 
will  rapidly  deteriorate  and  customers  will  become  disillusioned.  Allied  to 
this,  is  the  need  for  the  institution  to  be  able  to  develop  a  realistic  policy 
for  connection  of  consumers  to  the  system,  and  a  charging  policy,  in  or¬ 
der  to  enhance  revenue,  serve  the  maximum  number  of  people,  and 
maintain  affordability...  There  is  little  question  that  the  majority  of  peo¬ 
ple  are  willing  to  pay  for  safe,  reliable  drinking  water,  however  their  abil¬ 
ity  to  do  so  will  depend  on  macro-economic  conditions  unrelated  to  the 
project.  [Another]  key  lesson  is  that  macro-economic  conditions  can  de¬ 
teriorate  rapidly  through  earlier  unforeseen  circumstances,  jeopardizing 
the  sustainability  of  a  project. 

In  addition,  a  flood  that  occurred  in  1994  destroyed  some  of  the  infrastruc¬ 
ture,  handicapping  the  ability  of  the  utility  to  provide  (and  demand  pay¬ 
ment  for)  consistent  services  during  that  year. 

According  to  the  2010  census  data  (Minnesota  Population  Center  2015), 
which  included  only  very  limited  sampling  of  urban  areas  in  this  part  of 
the  country,  no  urban  residents  in  the  state  of  Borno  had  piped  water  and 
100%  of  the  water  supply  was  provided  from  surface  water.  This  might  in¬ 
dicate  that  even  where  the  water  is  supplied  by  a  system  of  tanks  and  pipes 
from  a  water  treatment  plant,  it  is  distributed  at  the  neighborhood  level 
via  trucks  or  hand  carrying.  In  recent  years,  businesses  that  pump,  treat, 
and  package  groundwater  for  household  use  have  become  commonplace. 
This  water,  also  known  as  “sachet  water”  (Figure  4)  has  become  an  im¬ 
portant  resource  for  urban  household  use  in  Nigeria.  (Dada  2009). 
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Figure  4.  Example  package  of  sachet  water 
(https://udeozochibuzo.wordpress.com/2014/01/14/cash-spiner- 
polyethylene-or-sachets-printing-business/). 


2. 4.3. 2  Conceptual  model 

Despite  the  negative  World  Bank  assessment  of  the  state  of  water  treat¬ 
ment  facilities,  the  UTP  indicates  the  water  treatment  is  at  least  marginally 
functional.  It  also  identifies  a  water  treatment  plant  near  the  airport,  four 
water  towers,  and  six  reservoirs.  Of  the  UTP  tanks  of  undefined  function, 

12  tanks  (or  sets  of  tanks)  were  selected  for  the  model  as  proxy  water  tanks 
based  on  one  or  more  of  the  following  criteria:  nearby  building  function, 
location  in  higher  elevations,  proximity  to  a  reservoir,  and  distance  from 
main  roads.  There  is  one  known  pumphouse  near  the  airport  that  presum¬ 
ably  is  associated  with  the  water  treatment  plant  and  fills  the  nearby  tanks 
from  the  reservoir.  The  facility  was  presumed  to  treat  water  for  drinking 
and  not  wastewater,  given  that  the  census  indicates  that  latrines  are  the 
norm  for  urban  areas  in  Borno  state.  It  was  also  presumed  that  any  reser¬ 
voirs  hold  untreated  water  and  thus  are  not  represented  in  the  urban  water 
infrastructure  sector  modeling. 

The  conceptual  model  has  been  populated  with  the  two  water  treatment 
plants  which  rely  on  electricity  from  the  substations.  These  treatment 
plants  feed  a  set  of  four  municipal  fresh  water  storage  towers  via  a  strong, 
presumably  piped  link;  they  also  feed  12  proxy  municipal  water  storage 
tanks  via  a  weak,  presumably  truck-delivered  link.  These  tanks  and  towers 
in  turn  supply  schools,  medical  facilities,  residential  areas,  and  commer¬ 
cial  areas  with  water  via  weak  links.  To  represent  the  “sachet  water”  com¬ 
ponent,  a  set  of  10  proxy  vendor-operated  well  water  locations  was  located 
at  a  subset  of  commercial  locations  near  surface  water  or  otherwise  in  lo¬ 
cations  that  might  be  amenable  to  groundwater  pumping.  These  busi¬ 
nesses  rely  on  electricity  provided  by  their  own  or  privately  owned  diesel 
generators. 
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2.4.4  Food  sector 

2.4.4. 1  Background 

The  state  of  Borno  has  traditionally  been  a  gateway  for  trade  in  agricul¬ 
tural  products,  livestock,  and  textiles  between  northern  Nigeria  and  neigh¬ 
boring  countries  including  Cameroon,  Chad  Niger,  and  the  Central  African 
Republic.  In  recent  years,  Boko  Haram  fighters  in  the  area  have  blocked 
highways,  destroyed  infrastructure,  and  killed  drivers.  This  mayhem  is  se¬ 
verely  dampening  commerce  and  preventing  agriculturalists  living  in  vil¬ 
lages  and  cities  from  accessing  their  lands  to  plant.  In  March  2014,  there 
were  an  estimated  250,000  internally  displaced  persons  in  the  three 
northeastern  states  of  the  country  including  Borno,  largely  due  to  the  ac¬ 
tivities  of  the  insurgent  group  (Ahmed  and  Eckel  2014). 

Many  urban  residents  of  Maiduguri  own  farms  outside  the  urban  area  and 
commute  from  the  city  to  the  farms  for  planting  and  harvesting  and  then 
bringing  the  harvests  back  to  the  city  for  sale  in  markets.  Livestock  from 
northern  pastoral  areas  of  the  country  also  are  delivered  to  the  city,  which 
has  a  cluster  of  animal  markets  and  processing  facilities  along  the  north¬ 
eastern  route  into  the  city.  Urban  markets  have  been  a  major  target  of  in¬ 
surgent  suicide  attacks  in  recent  years. 

Given  the  irrigated  areas  outside  of  the  city  (Figure  5)  as  well  as  the  preva¬ 
lence  of  boreholes  in  rural  areas,  the  water  sector  is  known  to  be  intimately 
tied  to  the  agricultural  sector  of  the  economy.  However,  this  analysis  fo¬ 
cuses  only  on  water  provision  for  household  and  commercial  use  within 
the  urban  area  and  does  not  explicitly  link  the  water  sector  to  the  food  sec¬ 
tor  in  the  infrastructure  model. 
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Figure  5.  Imagery  shows  the  agricultural  area  northeast  of  Maiduguri's  urban  area 
during  April  2014  (top)  and  November  2013  (bottom).  Fields  and  orchards  that  are 
green  during  the  dry  season  have  year-round  access  to  irrigation, 
either  via  wells  or  via  surface  water. 


2. 4. 4.2  Conceptual  model 

In  the  current  conceptual  model,  the  food  sector  is  only  minimally  repre¬ 
sented  with  three  livestock  processing  facilities  that  are  weakly  linked  to 
eight  major  market  locations  across  the  city.  In  turn,  these  market  loca¬ 
tions  support  residential  and  commercial  community  users  via  weak  links. 

2.4.5  Community  sector 

2.4.5. 1  Background 

Population  estimates  for  the  city  of  Maiduguri  vary  widely;  a  2011  estimate 
based  on  2006  census  figures  and  enumeration  areas  demarcated  by  the 
National  Population  Commission  bracketed  the  population  of  the  com¬ 
bined  Maiduguri  and  Jere  local  government  areas  between  approximately 
1  and  2.7  million  people  (Waziri  2009).  Since  that  enumeration,  however, 
an  estimated  250,000  Nigerians  have  been  forced  to  flee  their  homes  in 
northeastern  states  in  response  to  threats  and  attacks  associated  with 
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Boko  Haram.  Many  of  these  IDPs  have  fled  to  major  urban  areas  like  Mai- 
duguri,  so  the  current  population  is  likely  to  be  much  higher  than  esti¬ 
mates  made  in  the  previous  decade. 

Higher  population  density  occurs  in  several  districts:  Gwange,  Shehuri 
North,  Old  Maiduguri,  and  Yerwa,  each  of  which  has  over  10,000  people 
per  square  kilometer.  The  Ministry  of  Land  and  Survey  has  categorized 
housing  density  into  high-,  medium-,  and  low-density  areas  based  on  plot 
sizes,  although  this  information  is  not  always  well  correlated  with  the 
number  of  households  living  on  a  particular  plot.  The  number  of  houses  in 
a  given  administrative  ward  has  been  counted  in  a  2011  polio  eradication 
program,  and  independent  researchers  have  also  counted  individual 
houses  per  unit  area  in  select  wards  of  the  city  (Waziri  2009). 

The  northern  regions  in  Nigeria  are  populated  by  mainly  Muslim  ethnic 
groups  such  as  the  Hausa,  Fulani,  and  Kanuri,  whereas  populations  in  the 
south  of  the  country  are  largely  Christian  Yoruba  and  Igbo  (Ahmed  and 
Eckel  2014).  Maiduguri  is  majority  Muslim,  but  as  a  state  capital  and  uni¬ 
versity  town  also  is  has  some  Christian  presence;  the  UTP  Key  Features 
layer  identifies  over  90  mosques  and  fewer  than  10  churches.  In  the  model 
prepared  for  this  work,  however,  there  is  no  differentiation  between  com¬ 
munities  based  on  any  presumed  ethnic  or  religious  identities. 

2.4. 5. 2  Conceptual  model 

In  the  model,  rather  than  compiling  information  from  the  various  sources 
listed  above,  the  BTZ  layer  of  the  UTP  was  used.  These  BTZs  delineate  sets 
of  polygons  associated  with  different  land  use  and  residential  income  char¬ 
acteristics.  There  are  92  middle-income  residential,  30  poor  residential,  2 
shantytown  residential,  137  commercial,  5  industrial,  and  42  institutional 
polygons.  Together  these  polygons  cover  a  substantial  majority  of  the  ur¬ 
ban  area. 

The  population  of  each  community  polygon  was  estimated  by  using  the 
LandScan™  2010  dataset  created  by  Oak  Ridge  National  Laboratory 
(ORNL  2010).  To  produce  the  estimate,  the  LandScan  raster  was  con¬ 
verted  to  a  polygon  layer,  then  intersected  it  with  each  of  the  BTZ  layers, 
and  multiplied  it  by  the  fractional  area  of  each  community  within  that 
LandScan  grid  cell.  By  summing  this  value  over  all  LandScan  cells  in  a 
given  community  polygon,  an  estimate  was  obtained  for  the  populations 
represented  by  each  community.  These  community  population  values  can 
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be  used  in  the  network  equilibrium  model  for  calculations  of  demand  and 
queuing  time.  As  an  alternative,  another  simpler  representation  of  popula¬ 
tion  can  be  used  whereby  fixed  numbers  are  assigned  to  communities  by 
community  type,  with  the  mixed  commercial-residential  areas  found  near 
the  center  of  the  city  having  the  highest  populations  and  areas  more  com¬ 
monly  found  further  from  the  city  center  are  assigned  lower  values. 

The  BTZ  Institutional  community  layer  was  not  used  directly  in  this  analy¬ 
sis.  Institutions  in  this  layer  included  police,  court,  and  military  facilities 
as  well  as  universities,  but  the  type  of  access  cost  changes  modeled  with 
network  equilibrium  techniques  does  not  transfer  easily  to  the  types  of  re¬ 
sources  and  services  supplied  by  these  institutions;  a  community  near  a 
military  base  is  not  necessarily  better  positioned  to  have  better  “access”  to 
military- derived  resources.  Institutions  like  the  military  and  the  police 
which  “go  out”  into  the  communities  in  very  directed  ways  could  poten¬ 
tially  be  modeled  in  this  system,  but  would  require  a  third  type  of  linkage 
that  is  different  from  both  the  strong  and  weak  supports  currently  being 
used  for  other  resource  types.  This  third  type  of  linkage  would  itself  re¬ 
quire  very  specific  information  as  to  the  communities  targeted  for  specific 
operations,  as  well  as  some  way  to  quantify  the  value  of  the  provided  re¬ 
sources.  This  type  of  model  development  is  beyond  the  scope  of  the  cur¬ 
rent  project,  however,  which  aims  to  assess  the  effects  of  physical 
infrastructure  changes  on  local  populations. 

Therefore,  the  representative  institutions  chosen  were  the  84  schools  and 
4  hospitals  or  medical  facilities  identified  in  the  Key  Features  layer,  which 
were  then  treated  as  separate  infrastructure  components  that  served  the 
residential  populations.  Currently,  this  model  does  not  distinguish  be¬ 
tween  types  of  hospitals  or  levels  of  schools,  but  such  a  distinction  could 
be  included  in  future  analyses.  Of  primary  importance  for  the  community 
layers  is  defining  from  which  infrastructure  resources  they  pull  via  weak 
links  to  the  system.  In  this  case  study,  all  types  of  communities  needed  to 
access  water,  food,  education,  and  medical  resources,  with  the  exception  of 
the  few  industrial  locations  for  which  only  water  demand  was  represented 
(in  the  absence  of  more  specific  information  on  the  industrial  needs  at 
these  locations). 
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2.4.6  Transportation  sector 

2.4.6. 1  Background 

Residents  of  and  travelers  to  Maiduguri  use  a  wide  variety  of  transport  op¬ 
tions  from  cars,  trucks,  and  buses  to  motorbikes  and  bicycles,  to  donkeys 
and  camels,  and  to  pedestrian  transport.  The  city  also  has  an  airport  with 
flights  connecting  to  a  few  major  Nigerian  cities.  Motorbikes  in  particular 
have  become  common  and  in  some  cases,  they  dominant  in  congested  ur¬ 
ban  areas  in  Nigeria,  particularly  for  lower-  and  middle-class  residents. 

In  Maiduguri,  the  transportation  network  and  options  have  at  times  been 
significantly  affected  by  the  Boko  Haram  insurgency  and  the  military  re¬ 
sponses  to  it.  On  at  least  two  occasions,  the  government  instituted  a  tem¬ 
porary  ban  on  motorbikes,  apparently  with  varying  levels  of  success,  to 
address  violent  attacks  by  gunmen  from  that  particular  transportation 
platform  (“City-Wide  Motor  Bike  Ban”  2011).  On  certain  other  occasions 
such  as  Eid  al-Adha  observance  in  2015,  the  government  has  gone  so  far  as 
to  ban  all  cars,  public  transport,  horses,  donkeys,  and  camels  in  the  entire 
state  for  the  duration  of  the  holiday.  Vehicular  transit  in  and  out  of  the  city 
at  other  times  has  also  been  restricted  (“Nigeria’s  Boko  Haram  Crisis” 
2015).  These  events  highlight  the  notion  that  in  a  city  like  Maiduguri, 
where  the  population  relies  heavily  on  “weak  support”  via  the  transporta¬ 
tion  network  for  many  aspects  of  resource  acquisition,  it  may  be  important 
for  the  infrastructure  model  to  have  the  capability  to  represent  a  transpor¬ 
tation  network  that  is  dynamic  both  in  its  physical  structure  and  in  the 
cost  values  assigned  to  traversing  each  node. 

2.4.6.2  Conceptual  model 

To  represent  the  transportation  network  in  the  conceptual  model,  an 
ESRI  multipoint  dataset  including  over  9,500  road  segments  was  derived 
from  OSM  data.  Road  types  represented  in  the  dataset  included  primary, 
secondary,  tertiary,  residential,  and  trunk  segments.  Residential  segments 
were  removed  in  order  to  build  the  road  network  on  only  major  roads. 
Some  apparent  gaps  that  led  to  incorrect  dead-ends  in  the  resulting  net¬ 
work  were  cleaned  up  manually,  and  the  resulting  nodes  were  linked  to¬ 
gether  in  a  simplified  network  that  attached  each  node  to  the  four  closest 
nodes  within  a  threshold  distance  of  0.02  degrees,  thus  making  it  possible 


*  ESRI  develops  and  markets  software  that  melds  maps  with  data  for  expert  analysis. 
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to  represent  intersections.  This  simplified  road  network  was  developed 
primarily  with  a  goal  of  fast  processing  time;  in  a  full  model  run,  it  would 
be  desirable  to  use  a  full  road  network  with  appropriate  network  cost  val¬ 
ues  on  each  road  segment  according  to  the  different  road  types,  typical 
transit  speeds,  or  potentially  user-dependent  transit  speeds. 

2.5  Translation  of  conceptual  model  to  network  model  input 

After  the  conceptual  model  was  defined,  it  was  translated  (encoded)  into 
appropriate  input  for  the  network  equilibrium  analysis.  First,  the  relevant 
point  data  was  created  and  loaded  into  an  ESRI  geodatabase  with  feature 
datasets  representing  the  infrastructure  sectors  and  point  feature  classes 
representing  infrastructure  components  within  each  sector.  Then,  for  each 
component  node,  a  unique  node  identification  number  was  assigned  that 
encodes  the  sector  and  component  type  as  well  as  an  individual  compo¬ 
nent  number.  This  unique  and  succinct  identification  number  can  then  be 
used  both  by  humans  and  automated  systems  to  understand  and/or  pro¬ 
cess  the  network  information.  For  each  node,  an  attribute  was  created  to 
describe  operational  status  (a  binary  value  indicating  whether  or  not  the 
node  is  functional).  Next,  the  strong  and  weak  network  linkages  between 
nodes  were  defined  as  follows:  for  direct  support,  nodes  that  the  current 
node  directly  pushes  resources  to  in  the  form  of  strong  support  are  enu¬ 
merated;  for  weak  support,  nodes  that  the  current  node  can  receive  re¬ 
sources  from  are  enumerated.  An  additional  attribute  is  then  added  for  the 
transport  layer  to  represent  the  transport  nodes  directly  connected  to  the 
current  node.  In  this  process,  the  conceptual  model  is  used  as  a  reference 
to  develop  and  document  the  code  that  populates  the  geodatabase  with  the 
desired  attributes  that  define  the  network.  When  the  process  is  complete, 
information  on  all  nodes  in  the  geodatabase  is  exported  and  passed  to  the 
network  model.  This  dataset  includes  the  geographic  location  of  each 
point,  the  node  identification  numbers,  the  operational  status,  and  the 
linkage  information.  The  code  that  populates  the  geodatabase  serves, 
therefore,  as  the  final  documentation  for  the  network  structure. 

2.6  Alpha  model  development 

This  section  describes  the  analysis  framework  used  for  the  initial  Alpha 
layered  network  model,  addressing  mechanisms  of  community  behavior, 
mechanisms  of  disruption  propagation,  representations  of  system  equilib¬ 
rium,  and  calculations  of  social  impact  measures.  The  heuristic  method 
used  to  determine  system  equilibrium  is  also  explained  in  this  section. 
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2.6.1  Community  behavior 

Typically,  especially  after  disasters,  people  access  life-supporting  resources 
by  utilizing  the  transportation  layer.  The  trivial  case  that  people  get  re¬ 
sources  within  the  community,  (e.g.,  tap  water  at  each  household),  can  be 
incorporated  as  traveling  a  distance  of  zero.  Under  disruption,  fewer  re¬ 
source  nodes  remain  functioning,  which  can  lead  to  overwhelming  de¬ 
mand  concentration  and  associated  long  queuing,  longer  travel  detours, 
and  associated  traffic  congestion.  This  situation  significantly  increases 
people’s  resource  access  costs  and  may  further  lead  to  demand  loss  (i.e., 
people  giving  up  service)  if  the  costs  exceed  affordability  for  a  given  com¬ 
munity.  A  variant  of  the  traffic  equilibrium  model  was  implemented  to 
capture  this  issue. 

Suppose  the  demand  for  resource  r  for  each  community  i  e  Ic  is  inde¬ 
pendent  of  each  other.  Demand  is  denoted  by  d.  ,  and  the  corresponding 
resource  access  incurs  cost  c] .  Under  infrastructure  disruptions,  the  pro¬ 
curement  of  resources  normally  becomes  more  difficult  (e.g.,  obtaining 
water  from  a  distant  tank  rather  than  from  a  local  tap  or  pipeline  in  the 
case  of  a  water  network  disruption).  This  will  lead  to  the  increase  ofc'' , 

and  hence  the  decrease  of  d\ ,  that  is, 

d-  -  /t  (c;  u/;").  (i) 


where: 

D'  (-)  is  a  decreasing  inverse  demand  function  capturing  the 
elasticity 

d.°  is  the  initial  demand  without  infrastructure  failures. 

In  reality,  the  calculation  of  c.  is  complicated.  To  this  end,  the  following 
assumptions  are  considered: 

First,  each  community  i  can  access  resource  r  at  any  resource  node  j  e  /,. 
with  a  resource  access  cost 


Cij  '  Cij,i  Ci],s  ’ 


(2) 
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where: 

cij  t  is  the  travel  time  along  the  shortest  path  from  i  to  j  via  the 

transportation  layer 
cij  s  captures  service  (queuing)  time  at  j 

p  e  Py  is  a  path  (i.e.,  jDc/J  connecting  i  and  j  through  the 

transportation  layer  and 

gh,ke  It  be  the  demand  flow  on  transportation  node  k . 

Then,  the  result  is 


p£pi. 


(3) 


kep 


where: 


Ugk)  =  Sk 


1  +  a, 


rgpf‘ 


is  the  conventional  Bureau  of  Public 


Roads  (BPR)  function  (Sheffi  1985),  and 
ak  >  Pk  >  Yk  j  $k  are  node-specific  parameters. 

Furthermore,  the  resource  procurement  at  node  j  is  considered  to  line  up 
as  an  M/M/i  queue  (other  options  are  discussed  in  An  et  al.  2015).  Hence, 
the  average  waiting  time  is  computed  as: 


1 


(4) 


where 

A j  is  the  given  service  capacity  and 

hj  <  A j  is  the  total  acceptable  demand  flow  to  j . 

Second,  community  i  should  choose  the  resource  node(s)  with  the  mini¬ 
mum  access  cost,  i.e.: 


c.  =  mine... 

1  j(=j  V 


(5) 
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Given  ff  as  the  optimal  demand  allocation  along  path  p  produces  the  fol¬ 
lowing: 


X  (6) 

jsIr,psPy 

Z  ff  =  hj  ’  v/  G  f  ,VreR,  (7) 

and 

E  ff=g^kal,.  (8) 

Here,  note  that  the  scarcity  of  each  resource  is  modeled  by  its  access  diffi¬ 
culty  rather  than  by  the  amount  procured.  That  is,  the  supply  at  each 
working  resource  node  is  assumed  to  be  unlimited. 

Given  the  seven  infrastructure  layers,  the  transportation  layer  reaches  the 
above  equilibrium  (i.e.,  equations  (i)-(8))  when  no  community  can  reduce 
its  resource  access  cost  by  unilaterally  changing  its  resource  nodes  or  the 
travel  path  choices. 

2.6.2  Disruption  propagation  mechanism 

Now,  consider  the  scenario  where  there  are  disrupted  infrastructures. 

First,  to  maintain  working  status,  some  infrastructure  nodes  need  resource 
supply,  which  are  also  transported  via  the  transportation  layer  (e.g.,  a  die¬ 
sel  generator  accesses  fuel  from  diesel  tanks  filled  by  tanker  trucks). 
Therefore,  the  traffic  congestion  under  the  community  behavior  may  sig¬ 
nificantly  delay  or  block  resource  procurement  as  well,  or  even  cause  fail¬ 
ure  to  the  infrastructure.  This  type  of  failure  is  referred  to  as  a  resource 
failure. 

Second,  infrastructure  disruption  could  also  be  due  to  direct  and  obvious 
dependencies  based  on  physical  connections.  For  example,  a  water  tank 
with  an  electrical  pump  must  be  supported  by  a  nearby  power  source  (e.g., 
a  power  grid).  Such  a  dependency  is  normally  established  for  the  long  run, 
at  a  large  setup  cost,  and  hence  difficult  to  modify.  This  type  of  failure, 
caused  by  supporting  infrastructure  disruption,  is  referred  to  as  support 
failure. 
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Comparing  the  two  types  of  failure,  support  failure  is  found  to  be  one  ma¬ 
jor  cause  of  cascading  disruptions.  In  particular,  a  network  with  a  very 
dense  presence  of  links  (e.g.,  a  tree)  can  be  very  vulnerable— the  failure  of 
the  root  node  will  propagate  (cascade)  and  disrupt  the  entire  system.  In 
contrast,  resource  failure  normally  poses  a  potential  risk,  but  one  that  is 
less  obvious  to  foresee  in  advance.  However,  resource  failure  can  also 
bring  devastating  damage  to  the  entire  system,  especially  when  it  happens 
at  some  critical  node  that  can  cause  consequential  cascading  support  fail¬ 
ures. 

To  model  infrastructure  cascading  disruption,  let  x;  e  {0,1}  be  the  working 
status  of  each  node  i ,  where  xt  =  1  indicates  the  node  is  working,  or  o  indi¬ 
cates  otherwise.  Resource  failure  is  formulated  first.  Suppose  the  notations 
and  equilibrium  conditions  for  community  behavior  also  hold  for  a  generic 
node  iel\lc.  Let  e.  be  the  minimum  requirement  of  resource  r  to  run 

the  infrastructure  node  i .  Then,  if  there  exists  any  reR  such  that  d\  <  e' , 
then  a  resource  failure  happens  at  i .  As  for  the  support  failure,  Sj  c  /  is 
defined  to  be  the  subset  of  nodes  that  support  i .  That  is,  there  is  physical 
interdependency  between  each  node  in  Sj  and  i  such  that  if  any  i'  e  Sf 
fails,  it  would  lead  to  the  disruption  of  i . 

Therefore,  the  condition  for  node  i  to  be  functioning  can  be  expressed  as 
the  following: 


A=n*,nw  (9) 

jeSt  reR 

where  1M  is  the  indicator  function,  taking  the  value  of  l  if  the  condition  (•) 
holds  and  o,  otherwise.  Note  that  when  e.  -  0  holds,  it  implies  node  i  does 
not  need  resource  r .  For  notation  consistency,  for  community  i  e  /  . ,  set 
e.  =  0  for  all  r  and  Sj  =  0 ,  without  loss  of  generality.  This  notation  will 
guarantee  xt  =  1  always  holds  for  a  community  node;  i.e.,  no  community  is 
allowed  to  disappear. 

2.6.3  System  equilibrium 

With  the  resource  equilibrium  conditions  and  the  infrastructure  disrup¬ 
tion  mechanism,  the  system  equilibrium  is  now  ready  to  show.  Given  the 
collection  of  all  system  parameters  as 
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e:={a„/?„r„5»,/l „<,D, 


keL.iel  ,reR 


then 


T (0)  :=  jv; ,  d. , <  :  s.t.  Constraints  ( 1)-(9)|  , 

defines  all  system  variables  that  form  an  equilibrium  state.  Hence  s  e 
depicts  one  of  such  equilibrium  state.  Particularly,  let 

s  .-  {s  e  r(0) .  xt  - 1,  V/  e  If  SyStem  equilibrium  in  the  normal  sce¬ 

nario. 


It  can  be  observed  that  Constraints  (i)-(8)  guarantee  an  unique  equilib¬ 
rium  for  any  given  {  xf };  however,  it  is  possible  that  there  are  multiple 

equilibrium  states  that  have  different  working  infrastructures  {  xf }  and  yet 

satisfy  Constraints  (l)— (9).  Each  infrastructure  is  assumed  to  remain  func¬ 
tioning  unless  failure  is  absolutely  necessary,  and  hence  only  the  equilib¬ 
rium  with  the  most  working  infrastructures  is  considered.  Then,  given 
initially  disrupted  infrastructures,  (_)c/  such  that  xt  =  0,  i  e  Q ,  the  system 

equilibrium  after  failure  propagation  is  determined  by 


sy  :=  argmax 

ss{r(0):x,=o,  Vieg}  ig/ 


(10) 


2.6.4  Social  impact 

This  section  quantitatively  evaluates  the  social  impact  of  system  failure.  In 
general,  the  social  impact  is  multi-faceted.  First,  from  a  societal  perspec¬ 
tive,  the  reduction  of  demand  in  the  served  community  may  reflect  the  de¬ 
gree  of  people’s  unmet  needs,  suffering,  or  risk  of  dying  under  disruption 
situations.  Second,  from  an  economic  perspective,  the  resource  procure¬ 
ment  cost  significantly  increases;  meanwhile,  the  disrupted  infrastructure 
may  require  significant  economic  resources  and  effort  in  order  to  recover 
functionality.  To  this  end,  the  following  three  measures  are  considered  as 
proxies  for  evaluating  social  impact  across  a  set  of  communities  for  each 
type  of  resource:  (1)  total  demand  loss,  (2)  average  procurement  cost  in¬ 
crease,  and  (3)  total  infrastructure  failure. 

Given  the  initial  infrastructure  disruption  Q ,  the  system  equilibrium 
moves  from  the  normal  scenario  s°  to  se .  Hence,  for  each  resource  r  e  R , 
the  total  demand  loss  is: 
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ielr 


the  average  procurement  cost  increase: 


and  the  total  infrastructure  failure: 


A*=Z(4-44 


(id 


(12) 


(13) 


More  sophisticated  impact  measures  are  possible,  but  they  are  beyond  the 
scope  of  this  paper  and  are  left  for  future  research. 

2.7  Solution  approach 

In  this  section,  a  heuristic  algorithm  is  developed  that  extends  the  conven¬ 
tional  Frank-Wolfe  algorithm  (Sheffi  1985)  to  obtain  se .  Noticing  the  re¬ 
source  demand  is  elastic  and  there  exists  a  queuing  problem  at  the 
resource  node,  an  equivalent  problem  with  constant  demand  in  a  new  net¬ 
work  representation  was  constructed. 

In  the  new  network,  a  virtual  node  f  is  added  for  each  resource  r  eR  . 
Then,  two  types  of  virtual  roads  (nodes)  are  added  to  the  transportation 
layer  to  capture  demand  elasticity  and  queuing  issues,  respectively.  First,  a 
virtual  road  kj  is  added  from  any  j  e  /.  to  f ,  with 


TAgk) 


(14) 


and  the  original  demand  flow  to  j  is  directed  into  a  new  destination  ir . 
Hence,  gk  exactly  takes  the  value  of  h  j ,  and  the  queuing  cost  is  converted 

to  an  equivalent  transportation  cost.  When  Xj  =  0 ,  the  transportation  cost 
goes  to  infinity,  which  guarantees  no  demand  flow  goes  to  j .  Second,  for 
each  i'  ,r  eR,  a  virtual  road  k ■  is  added  from  any  i  e  /  to  ir ,  with 
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T,r  («„)  =  ■ 4'  («„).  (15) 

where  D'  is  the  inverse  function  of  D"  as  in  (l),  and  fix  the  resource  de¬ 
mand  at  i  to  be  d.° . 

In  so  doing,  all  resource  demands  at  i  are  constant  regardless  of  the  re¬ 
source  procurement  cost  or  working  status  of  infrastructures  (overflows  or 
lost  demand  are  “sent”  via  virtual  links).  Furthermore,  all  demand  of  re¬ 
source  r  goes  to  the  same  node  f .  As  such,  the  problem  has  been  con¬ 
verted  to  a  conventional  network  equilibrium  problem. 

Then,  the  following  iterative  algorithm  is  proposed  to  estimate  se : 

Step  o.  Initialization.  Regardless  of  resource  failure,  iteratively  add  i  to  Q 
if  Sj  f)  Q  ^  0 ,  until  no  new  nodes  can  be  added;  then,  set  x([0]  =  0  for  all 
i  g  Q  and  l,  otherwise;  next,  conduct  traffic  assignment  based  on  x,.[0]  to 
obtain  gA  [0]  for  all  k  e  It ; 

Step  l.  Ti'affic  assignment.  Conduct  traffic  assignment  based  on  xfn]  to 
obtain  demand  flow  g[ ,  and  update  gk [n  + 1]  =  (1  -  pn )gk [n]  +  png[ ; 

Step  2.  Update  resource  failure.  Construct  Q’  to  contain  all  nodes  i , 
where  there  exists  r  e  R ,  such  that  d'f  -  g,  [n  + 1]  <  e. ; 

Step  3 .  Renew  failure  propagation.  Iteratively  set  xfn  + 1]  =  0  for  i  e  Q'  or 
.  D  Q'  ^  0 ,  and  expand  Q'  =  Q'U  {/} ,  until  no  nodes  can  be  added; 

Step  4 .  Stopping  rule.  If  n  =  N  or  ^  |  gk [//  +  l]-g/(  [n]  \<d ,  terminate,  oth- 

k<=It 

erwise,  set  n  =  n  + 1  and  repeat  Steps  l  to  3. 

00 

Here,  pn  is  a  decreasing  search  step  such  that  lim  pn  =  0  and  V  pn  -o-  +co , 

n= 0 

N  is  the  maximum  iteration  number,  and  £  is  the  small  positive  thresh¬ 
old. 
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2.8  Alternative  model  formulations 

In  addition  to  the  original  model  formulation,  there  are  two  additional 
ways  in  which  the  model  has  been  modified  in  order  to  address  specific 
end-user  questions: 

1.  How  can  the  model  be  modified  to  reflect  certain  aspects  of  uncer¬ 
tainty,  such  as  uncertainty  in  exactly  how  fragile  a  given  node  is  under 
disruption? 

2.  After  some  infrastructure  nodes  have  been  damaged,  how  can  this  net¬ 
work  be  used  to  decide  which  nodes  need  to  be  restored  to  maximize 
social  impact? 

These  questions  are  addressed  via  specific  modeling  techniques  as  ex¬ 
plained  in  the  subsections  that  follow. 

2.8.1  Probabilistic  failure  using  Monte-Carlo  simulations 

An  end-user  might  want  to  know  how  to  represent  a  case  where  there  is  in¬ 
complete  information  about  the  true  fragility  of  infrastructure  under  the 
given  failure  propagation  mechanisms.  In  this  modification,  rather  than 
presuming  that  infrastructure  failure  occurs  completely  and  consistently 
when  certain  predefined  failure  conditions  are  met,  instead  probabilistic 
failure  analysis  is  used,  whereby  each  failure  propagation  occurs  based  on 
a  certain  predefined  probability  of  failure.  This  technique  can  help  model 
results  under  incomplete  information  because  the  mean,  range,  and  distri¬ 
bution  of  the  results  of  multiple  model  runs  can  be  analyzed  to  understand 
the  implications  of  uncertainties  in  input  failure  probabilities. 

In  the  Monte-Carlo  case  study,  a  failure  probability  for  all  facilities  of  0.9 
was  used,  and  300  model  iterations  were  run.  In  a  single  iteration,  all  fa¬ 
cilities  have  status  values  that  are  either  o  or  1  for  representing  nonfunc¬ 
tioning  or  functioning.  However,  when  these  iterations  were  combined 
and  the  mean  status  was  used  to  represent  overall  model  behavior,  it  was 
noted  that  what  counts  as  failure  needed  to  be  defined,  given  that  mean 
status  results  would  no  longer  be  represented  as  integer  values  of  o  and  1. 
The  number  of  failed  facilities  in  each  sector  were  counted  under  the 
Monte-Carlo  simulation  as  if  a  final  mean  status  constitutes  failure  if  it  is 
less  than  0.5. 
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2.8.2  Restore,  recover,  and  rebuild  using  the  genetic  algorithm 

In  the  “Restore-Recover- Rebuild”  case  study,  the  team  was  given  an  al¬ 
ready-disrupted  system  and  a  set  of  nodes  that  could  be  restored  or  even 
rebuilt  in  new  locations.  The  infrastructure  network  model  was  viewed  as  a 
function  that  maps  each  subset  of  possible  combinations  of  nodes  to  be  re¬ 
stored  to  a  set  of  social  welfare  indicators,  such  as  the  cumulative  time  for 
all  communities  to  access  water,  markets,  healthcare,  etc.  A  method  can 
then  be  set  up  to  explore  the  best  possible  combination  of  nodes  to  restore, 
given  that  aim  is  to  maximize  social  benefits  while  only  restoring/rebuild¬ 
ing  a  certain  fixed  number  of  nodes  at  any  point  in  time.  This  method 
translates  mathematically  to  a  nonlinear  non-convex  optimization  prob¬ 
lem,  which  can  be  solved  approximately  but  efficiently  by  using  a  genetic 
algorithm.  In  this  method,  rather  than  performing  an  exhaustive  and  com¬ 
putationally  expensive  calculation  that  runs  all  possible  restoration  sce¬ 
narios  and  compares  all  resulting  social  welfare  indicators  to  find  the 
optimum  solution,  the  following  steps  were  taken: 

•  select  a  smaller  population  of  candidate  solutions  (node  combinations 
to  restore), 

•  run  the  model  on  them  and  select  the  best-performing  solutions, 

•  modify  the  solutions  by  crossover  (random  mixing  of  parts  of  two  good 
solutions)  and  mutation  (random  perturbation  of  a  solution), 

•  run  the  model  on  those  solutions  again,  and 

•  use  the  resulting  social  welfare  indicators  to  discard  underperforming 
modifications,  returning  the  more  successful  modifications  to  the  pop¬ 
ulation  pool. 

The  process  is  then  repeated  a  certain  number  of  times  to  explore  the 
problem  space.  The  process  can  be  modified  to  optimize  solutions  for  a 
particular  social  indicator,  a  particular  neighborhood,  or  an  area  of  the  city 
most  of  interest  to  an  end  user. 

In  the  results  below,  the  case  is  shown  where,  after  the  disruption  of  all 
medical  facilities,  the  effects  of  restoration  of  an  existing  medical  facility 
were  compared  with  the  construction  of  new  facilities.  The  pool  of  poten¬ 
tial  locations  for  new  construction  was  constrained  to  the  locations  of  ex¬ 
isting  community  nodes. 
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2.9  Beta  model 

A  number  of  other  modifications  were  now  made  to  the  original  (Alpha) 
model  to  account  for  common  behaviors  of  physical  and  economic  systems 
not  previously  represented.  The  second  (Beta)  model  version  that  ad¬ 
dresses  these  behaviors  is  described  below: 

•  Rather  than  a  binary  facility  status  of  fully  failed  or  fully  functional,  the 
facility  status  was  represented  along  a  continuum  from  o  to  l,  so  that 
minor  malfunctions  and  associated  reduced  capacity  can  be  repre¬ 
sented.  Let  x‘  e  be  the  working  status  of  each  node  z ,  where  the 

x  —  1 

node  is  fully  functioning  if  '  ,  completely  failed  if  xt  =  0,  and  par¬ 

tially  working  otherwise. 

•  A  “resource  limit”  or  “service  capacity”  was  defined,  such  that  no  facil¬ 
ity  can  be  presumed  to  offer  infinite  resources  during  a  given  model  it¬ 
eration.  Define  ^  >  ^  and  ~j  ~  ^  to  be,  respectively,  the  maximum  and 
minimum  service  capacity  of  facility  -1 . 

•  This  service  capacity  is  further  associated  with  facility  status,  such  that 
the  resource  limit  of  the  facility  scales  linearly  between  the  minimum 
and  maximum  capacity  values,  according  to  its  status: 

Aj  =  XjAj  +  (1  - Xj )Xj ,  V/  e  /,.,Vr  e  R  ( 16) 

•  In  the  Alpha  model,  although  queuing  cost  was  mathematically  repre¬ 
sented  as  an  M/M/ 1  queue,  the  fact  that  service  capacity  was  unlimited 
meant  that  this  queuing  cost  was  in  practice  always  zero,  and  the  total 
calculated  cost  was  a  function  of  transportation  costs  only.  Under  the 
new  consideration  of  limited  resources,  this  “queuing”  (which  can  rep¬ 
resent  the  rate  of  resource  production  to  resource  consumption  as  well 
as  the  idea  of  “waiting  in  line  to  be  served”)  can  now  correctly  result  in 
access  costs  associated  with  transportation  costs  and  queuing  costs. 

•  Furthermore,  the  elasticity  of  demand  was  represented  in  the  Beta 
model.  The  Alpha  model  presumed  that  demand  for  a  resource  was 
constant  below  a  certain  cost  threshold,  at  which  point  the  consumer 
could  no  longer  afford  the  resource,  and  demand  collapsed  completely. 
In  the  Beta  model,  demand  decreases  as  a  function  of  increasing  cost  to 
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represent  the  fact  that  people  often  make  do  with  less  as  costs  rise,  ra¬ 
ther  than  waiting  to  hit  some  cost  ceiling  and  then  being  forced  to 
make  do  with  nothing. 

•  In  order  to  compute  system  equilibrium  on  a  network  that  includes 
these  new  components  of  queuing  cost  and  demand  loss,  an  expanded 
network  (Figure  6)  was  introduced  with  additional  nodes  that: 

o  funnel  all  resource  acquisitions  through  “queuing  nodes”  that  track 
queuing  costs, 

o  provide  a  “sink”  node  linked  to  each  facility  seeking  a  resource,  to 
which  the  value  of  community  demand  is  passed,  and 
o  create  a  “dummy”  node  for  each  community  that  bypasses  the 
roads,  the  resource  locations,  and  the  queuing  nodes  to  track  the 
fraction  of  a  community  that  no  longer  uses  the  service  when  de¬ 
mand  falls. 

Computing  the  Nash  equilibrium  on  this  augmented  network  means  that  a 
flow  pattern  is  found  in  it,  such  that  each  resource  user  finds  its  best  strat¬ 
egy,  i.e.,  choosing  resource  retrieval  path  and  resource  supply  location  (or 
giving  up  the  need  for  resource),  and  facility  working  status  (and  hence 
service  capacity)  agrees  with  their  weak  support  input  based  on  its  incom¬ 
ing  resource  flow. 


Figure  6.  Augmented  network  representation. 
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3  Results 


In  this  section,  different  versions  of  the  conceptual  model  are  imple¬ 
mented  to  analyze  the  social  impact  of  infrastructure  failure  based  on  the 
conceptual  model’s  combination  of  real-world  and  inferred  data  from  the 
city  of  Maiduguri,  Nigeria. 

3.1  Alpha  model  results 

3.1.1  Case  Study  AO 

The  residential  communities  in  the  urban  area  of  Maiduguri  are  catego¬ 
rized  into  four  groups  based  on  their  incomes;  these  categories  are  middle- 
income,  mixed  middle-income  and  commercial,  poor-income,  and  shanty¬ 
town.  Their  geographical  distributions  are  shown  in  Figure  7,  overlaid  on 
population  density.  The  population  distribution  is  observed  to  be  clustered 
near  the  city  center  (where  the  population  has  mostly  higher  income),  and 
it  gradually  diminishes  to  the  peripheral  (where  the  population  has  mostly 
lower  income).  The  population  of  each  community  polygon  was  estimated 
by  using  the  LandScan20io  dataset  (ORNL  2010).  The  resource  demand 
for  community  i  e  Ic  was  further  estimated  by: 


d 


r  0 


df,  if  r  e  {water,  food} 

<  0.4 df ,  if  r  €  {education,  healthcare} , 
0,  otherwise 


(16) 


where  df  is  the  population  of  community  i .  The  factor  0.4  reflects  the  re¬ 
source  demand  for  education,  and  healthcare  is  relatively  less  critical  than 
water  and  food.  In  addition,  an  inverse  demand  function  of  resource  r  for 
community  i  was  considered  as: 


0, 


if<<  max  {0.6c;r,  10<°} 
otherwise 


(17) 


where  c'°  is  the  resource  procurement  cost  under  the  normal  scenario, 
and  c\  is  the  maximal  resource  procurement  cost  among  the  community 
in  the  same  category  as  i .  This  is  to  say,  a  community  with  higher  c'  °  is 
more  robust  under  disruption  since  its  people  are  relatively  accustomed  to 


ERDC  TR-16-11 


39 


the  high-cost  (time)  life  style.  On  the  contrary,  a  community  with  low  c'° 

becomes  relatively  sensitive  to  the  cost  increase  under  disruption.  Here,  to 
avoid  unrealism,  c\  is  introduced  to  represent  the  minimum  affordability 

to  communities  with  extremely  low  c[° . 


Figure  7.  Maiduguri  community  and  population  distribution. 
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The  infrastructure  locations  are  presented  in  Figure  8,  where  each  point 
indicates  an  infrastructure  node,  and  the  type/layer  is  differentiated  by  the 
marker  shape.  The  structure  of  the  seven  infrastructure  layers  are  de¬ 
scribed  in  detail  in  a  previous  section,  but  are  summarized  here  as  follows: 
(i)  the  electricity  layer  supports  the  water  layer,  and  they  jointly  support 
education  and  healthcare  layers;  (ii)  the  transportation  layer  neither  is 
supported  nor  supports  other  layers,  but  purely  provides  mobility  to  the 
system;  (iii)  the  food  layer  only  provides  food  to  the  community,  where  the 
downstream  food  market  node  seeks  supply  from  the  upstream  food  pro¬ 
cessing  facilities;  and  (iv)  the  fuel  layer  only  provides  fuel  to  the  infrastruc¬ 
ture  nodes  such  as  those  in  the  electricity  layer.  The  resource  demand  of 
the  infrastructure  node  normally  does  not  generate  huge  traffic,  hence 
d.°  =  1  was  assigned  for  each  infrastructure  node  that  would  be  seeking  re¬ 
sources  (e.g.,  fuel  and  food);  the  inverse  demand  function  of  an  infrastruc¬ 
ture  i  is 
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if  c\  <  max {0.6c'  ,1.2c'0} 
otherwise 


(18) 


where  the  factor  1.2  is  smaller  than  that  of  a  community.  This  is  because 
the  flexibility  to  change  a  resource  supplier  for  the  infrastructure  is  nor¬ 
mally  far  less  than  that  of  a  resource  supplier.  It  should  also  be  noted  that 
some  infrastructure  node  within  a  resource  layer  might  not  directly  pro¬ 
vide  resources  but  support  other  downstream  resource  node,  which  in¬ 
cludes  water  treatment  plants,  fuel  central  depot,  and  food  processing 
facilities. 


Figure  8.  Maiduguri  infrastructure  distribution. 
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The  rest  of  the  parameters  are  set  as  follows:  ak  =  0. 15  min,  /3k  =  4, 
yk  =  10,000  vehicles  per  lane  per  day,  Sk  is  computed  based  on  actual  dis¬ 
tance  between  nodes  at  a  free-flow  speed  of  35  mph,  Xk  is  set  to  be  infinity 
to  indicate  the  resource  supply  is  sufficient,  e,'  =  d.°  for  iel\lc,  and  for 
running  the  algorithm,  N  =  50  and  pn=H n. 

A  scenario  was  considered  wherein  disruption  initiates  at  the  only  power 
substation  and  then  the  failure  propagation  and  social  impacts  were  inves¬ 
tigated.  The  electricity  substation  is  a  critical  node  that  supports  a  variety 
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of  other  infrastructure  nodes,  including  water  treatment  plants  and  elec¬ 
tricity  transformers.  After  implementing  this  algorithm,  the  system  con¬ 
verges  to  an  equilibrium,  and  the  infrastructure  status  is  shown  as  Figure 
9.  As  expected,  all  electricity  transformers  are  shut  down,  the  electricity 
network  is  disabled,  while  only  local  electricity  generators  can  work  based 
on  fuel,  providing  very  limited  (and  expensive)  power  supply  to  nearby 
communities.  Furthermore,  the  water  network  is  also  disrupted,  while 
only  local  water  vendors  can  pump  water  from  wells.  As  a  result,  educa¬ 
tional  institutions  and  healthcare  facilities  are  all  disrupted  due  to  the 
shortage  of  power  and  water.  In  addition,  a  fuel  tank  fails  because  the  cost 
to  get  fuel  supply  from  the  central  fuel  depot  becomes  too  high  due  to  local 
congestion. 

Figure  9.  Status  of  infrastructure  components  after  a  cascading  disruption. 
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The  social  impacts  are  estimated  as  follows.  In  total,  64%  of  the  water 
nodes,  0%  of  the  food  nodes,  11%  of  the  fuel  nodes,  and  100%  of  education 
and  healthcare  nodes  are  disrupted  (i.e.,  A"ater  =  18 ,  A*od  =  0 ,  A“  =  1 , 

A  education  _  g4  ^  anc[  ^healthcare  _  4)  Interestingly,  after  disruption,  the  average 

access  cost  for  food  even  decreases  from  0.75  to  0.54.  This  is  mainly  due  to 
the  loss  of  education  and  healthcare  needs  which  relieve  traffic  congestion 
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in  some  parts  of  the  city.  The  overall  average  access  cost  for  water  is  al¬ 
most  unchanged  for  the  same  reason.  A  detailed  look  of  the  spatial  distri¬ 
bution  of  community’s  water  access  cost  is  shown  in  Figure  to.  Compared 
to  that  in  the  normal  scenario,  the  cost  in  the  southwest  region,  which  was 
previously  very  high  due  to  existing  congestion,  significantly  decreases  due 
to  the  reduced  travel  demand.  In  contrary,  in  other  regions,  the  cost 
slightly  increases  as  original  water  suppliers  were  likely  to  be  disrupted. 

Figure  10.  Community  water  costs  before  (top)  and  after  (bottom)  disruption. 
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3.1.2  Case  Study  A1 

For  a  second  case  study  scenario  in  Maiduguri,  a  situation  was  introduced 
in  which  failure  is  caused  to  one  of  four  water  towers  supplied  by  clean  wa¬ 
ter  via  the  water  treatment  plants.  Both  before  and  after  the  failure,  the 
user  equilibrium  algorithm  calculates  an  access  cost  distribution  over  all 
the  communities  that  require  the  resource  such  that  after  the  equilibrium 
is  calculated,  no  individual  user  can  unilaterally  lower  his  or  her  access 
costs  by  choosing  to  obtain  the  resource  from  a  different  location.  These 
access  cost  results  are  provided  at  the  point  locations  of  the  communities 
at  which  they  are  calculated.  These  points  can  then  be  interpolated  to  cre¬ 
ate  a  cost  map  (Figure  n).  This  image  shows  that  costs  are  lowest  for  com¬ 
munities  nearer  to  the  water  towers,  while  a  band  of  higher-cost  access  is 
observed  to  run  north-south  through  the  western  part  of  the  city.  The  cost 
value  calculated  is  a  relative  number,  not  intended  to  represent  monetary 
units. 

When  the  northernmost  water  tower  fails,  the  user  equilibrium  is  recalcu¬ 
lated,  and  the  ratio  of  the  access  costs  before  and  after  the  failure  can  be 
displayed  (Figure  12).  This  mapping  clearly  delineates  the  areas  that  would 
be  most  immediately  affected  by  the  loss  of  this  particular  infrastructure 
component.  The  locations  of  residential  communities  of  varying  income 
levels  are  overlaid  on  this  cost  ratio  map.  Note  that  many  residential  com¬ 
munities  within  the  area  most  affected  by  the  loss  of  infrastructure  are 
middle-income  areas  and  only  two  are  poor  or  shantytown  areas.  The  ma¬ 
jority  of  the  highly  affected  areas  are  mixed-use  commercial/residential 
areas  (not  shown  in  Figure  12  for  readability  reasons;  mixed-use  commu¬ 
nities  are  mapped  in  Figure  7). 
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Figure  11.  Initial  water  access  costs  in  case  study  Al. 
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Figure  12.  Water  access  cost  ratio,  final,  case  study  A1 . 
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3.1.3  Alpha  model  sensitivity  analysis 

The  above  Case  Study  Ao  (outlined  in  3.1.1)  is  denoted  as  a  benchmark 
Scenario  o,  and  it  now  examines  several  representative  scenarios  under 
different  parameter  settings.  The  corresponding  societal  impacts  are  sum¬ 
marized  in  Table  1. 
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First,  in  Scenarios  l  and  2,  the  flexibility  for  an  infrastructure  to  change 
resource  supplier  is  perturbed.  For  example,  the  coefficient  of  c.°  in  Equa¬ 
tion  18.  Scenario  1  suggests  that,  with  a  lower  flexibility  to  change,  infra¬ 
structures  are  more  vulnerable  to  disruption  caused  by  resource  flow 
congestion/queuing.  Compared  with  the  benchmark  scenario,  more  infra¬ 
structure  (including  five  water  vendors  and  three  food  markets)  are  dis¬ 
rupted  due  to  resource  failure.  Consequently,  the  average  food  access  cost 
increases  by  95%,  and  now  0.9%  of  the  population  has  lost  its  food  supply. 
Further,  the  average  water  access  cost  dramatically  increases— by  more 
than  five  times  that  in  the  normal  scenario— since  only  five  water  nodes 
now  remain  working.  These  results  imply  that  the  system  can  be  very  sen¬ 
sitive  to  the  infrastructure’s  flexibility  for  changing  suppliers.  This  result  is 
due  mainly  because  of  the  system’s  positive  feedback  to  the  traffic  conges¬ 
tion  (i.e.,  stubborn-to-change  infrastructures  are  likely  to  have  resource 
failure  and  intensify  the  traffic  congestion,  which  in  turn  further  prevent 
the  change  of  suppliers).  On  the  other  hand,  the  system  is  more  robust  at  a 
higher  infrastructure  flexibility,  as  shown  in  Scenario  2.  The  only  noticea¬ 
ble  change  is  the  survival  of  the  fuel  tank  that  was  disrupted  in  the  bench¬ 
mark  scenario. 

Next,  in  Scenarios  3  and  4,  the  impact  of  having  a  smaller  road  capacity 
is  considered.  The  results  are  quite  dramatic.  Compared  with  the  bench¬ 
mark  scenario,  10%  of  capacity  reduction  can  cause  severe  food  and  water 
cost  increases  and  demand  losses.  This  comparison  confirms  how  resource 
failures  can  sweep  down  the  infrastructure  system.  However,  it  should  also 
be  noted  that  the  social  impacts  at  10%  and  25%  capacity  reduction  are 
relatively  close  (e.g.,  there  exists  a  diminishing  marginal  impact).  This  di¬ 
minishing  impact  may  occur  because  most  communities  have  lost  service 
even  at  10%  capacity  reduction. 


Table  1.  Sensitivity  analysis  results. 


Scenarios 

Food 

Water 

Access 

Cost 

Increase 

Demand 

Loss 

#Node 

Failed 
(Total  11) 

Access 

Cost 

Increase 

Demand 

Loss 

#Node 

Failed  (Total 
28) 

0:  Case  Study 

-36% 

0.0% 

0 

0.5% 

0.0% 

18 

i.ic;0 

1:  1 

95% 

0.9% 

3 

531% 

55.2% 

23 

1  'j  rO 

2:  1.3C, 

-37% 

0.0% 

0 

0.5% 

0.0% 

18 
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Food 

Water 

Scenarios 

Access 

Cost 

Increase 

Demand 

Loss 

#Node 

Failed 
(Total  11) 

Access 

Cost 

Increase 

Demand 

Loss 

#Node 

Failed  (Total 
28) 

3:  90 %8k 

536% 

78.1% 

7 

465% 

56.2% 

23 

4:  75%Sk 

601% 

80.0% 

7 

531% 

58.3% 

23 

5:  Initial  Water 

3% 

0.0% 

0 

6% 

0.0% 

7 

6:  Initial  Fuel 

25% 

0.0% 

0 

100% 

0.0% 

14 

7:  Initial  Fuel  &  Water; 

1.3<° 

72% 

1.9% 

1 

149% 

7.6% 

17 

Scenario  5  explores  a  different  initial  disruption  at  a  water  treatment 
plant.  Counter-intuitively,  although  water  treatment  plant  is  the  key  infra¬ 
structure  in  the  water  layer,  shutting  it  down  does  not  significantly  affect 
the  entire  system.  Although  three  water  towers  (due  to  support  failure) 
and  four  water  tanks  (due  to  resource  failure)  are  disrupted,  the  water  ac¬ 
cess  cost  only  slightly  increases  by  6.0%,  and  there  is  no  water  demand 
loss.  This  is  mainly  because  ten  local  water  vendors  can  still  independently 
provide  water  to  the  population  despite  the  water  network  disruption.  Fur¬ 
ther,  no  disruptions  are  found  in  the  food  layer,  although  the  increase  in 
traffic  congestion  (due  to  water  shipments)  still  causes  a  3%  increase  in 
the  average  food  access  cost. 

In  Scenario  6,  it  is  assumed  that  the  initial  disruption  occurs  at  the  cen¬ 
tral  fuel  depot,  which  directly  cuts  down  the  fuel  supply  to  the  entire  sys¬ 
tem.  As  a  result,  ten  local  water  vendors  are  all  disrupted  (due  to  resource 
failure)  and  four  water  tanks  are  further  disrupted  due  to  conges¬ 
tion/queuing  (due  to  increased  water  shipment  flow).  However,  the  com¬ 
munity  can  still  survive  well  without  fuel  supply,  whereas  only  access  costs 
for  food  and  water  bear  some  modest  increases.  Therefore,  the  fuel  layer 
does  not  play  a  critical  role  in  the  infrastructure  system  if  the  power  grid 
and  water  network  still  function  well. 

Finally,  in  Scenario  7,  the  central  fuel  depot  is  still  considered  as  the  ini¬ 
tial  disruption  node,  but  the  infrastructures’  flexibility  to  change  suppliers 
is  increased  to  1.3.  The  result  appears  to  be  quite  counter-intuitive  at  first 
glance.  Even  more  infrastructures  are  disrupted  in  this  scenario  than  in 
Scenario  6.  This  “paradoxical”  observation  is  a  reflection  on  the  interesting 
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interdependency  due  to  self-interested  competitions  among  different  in¬ 
frastructures.  When  a  certain  amount  of  infrastructures  shut  down,  as  in 
Scenario  6,  other  unrelated  infrastructure  will  survive  better  as  conges¬ 
tion/queuing  becomes  mitigated  and  resource  procurement  becomes  eas¬ 
ier.  On  the  contrary,  the  survival  of  a  smaller  number  of  infrastructures 
may  boost  traffic  congestion  and  even  disrupt  unrelated  infrastructures 
that  are  more  sensitive  to  resource  failure. 

3.2  Alternative  model  results 

3.2.1  Probabilistic  failure  with  Monte-Carlo  simulation 

The  results  of  using  a  probabilistic  failure  scheme  with  a  uniform  failure 
probability  of  0.9  shows  that  after  300  iterations,  the  mean  number  of 
functioning  nodes  does  indeed  change  significantly  across  most  sectors 
from  the  deterministic  benchmark  case  (Table  2).  However,  it  is  clear  that 
being  able  to  establish  even  the  direction  of  that  change  depends  on  how 
failure  is  ultimately  defined.  Therefore,  it  is  recommend  that  more  care  be 
taken  to  establish  how  resource  failure  is  defined  within  the  model,  which 
is  addressed  further  in  Section  3.3. 

Table  2.  Number  of  facility  failures,  by  sector,  under  probabilistic  failure  scenario. 


Failure  of  Facilities 

(Failure  Probability  =  0.9,  300  iterations) 


Level 

Total  number 

Benchmark  case 

(deterministic) 

Monte-Carlo  simulation 
(count  as  failure  if 
survival  rate  <  0.5) 

Monte-Carlo  simulation 
(count  as  failure  if 
survival  rate  <  0.1) 

Total 

544 

320 

428 

167 

Power 

408 

213 

309  117 

Fuel  9 

1 

7  0 

Food  11 

0 

1  0 

Water  28 

18 

23  2 

Education 

84 

84 

84  45 

Medical  care 

4 

4 

4  3 
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3.2.2  Restore,  recover,  and  rebuild  using  the  genetic  algorithm 

Figure  13  (a-d)  show  the  costs  associated  with  access  to  medical  facilities 
under  different  situations:  (a)  initial  state,  before  disruption;  (b)  after  dis¬ 
ruption,  where  costs  increase  significantly;  (c)  costs  after  one  of  the  medi¬ 
cal  facilities  is  restored;  and  (d)  costs  after  the  genetic  algorithm  has 
evaluated  the  optimum  placement  for  two  new  hospitals  (without  rebuild¬ 
ing  the  old  ones),  in  which  lower  overall  costs  are  seen  across  the  city. 
Note  that  in  these  figures,  the  scales  are  not  equivalent  in  order  to  show 
the  geospatial  differentiation  produced  by  the  network  model  equilibrium. 


Figure  13.  Medical  facilities  access  costs  under  Repair-Restore-Rebuild  scenario,  (a) 
prior  to  disruption  (low  costs),  (b)  after  disruption  (high  costs),  (c)  after  restoration  of 
one  previously  existing  facility  (moderate  costs),  (d)  after  rebuilding  in  two  new 
locations  (relatively  low  costs)  as  optimized  using  the  genetic  algorithm.  Note  that 
image  colors  are  scaled  differently  between  the  scenarios,  but  always  with  light  blue 
representing  the  lowest  cost  and  bright  pink  representing  the  highest  cost. 
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Medical  Facilities  Cost- After  Rebuild 
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3.3  Beta  model  results 

3.3.1  Case  Study  BO 

As  in  Case  Study  Ao,  a  scenario  was  considered  where  disruption  initiates 
at  the  main  power  substation  in  the  center  of  the  city,  and  the  failure  prop¬ 
agation  and  social  impacts  are  then  investigated.  After  implementing  the 
Alpha  model’s  algorithm,  the  system  converges  to  an  equilibrium,  and  the 
resulting  infrastructure  status  is  shown  in  Figure  14,  where  a  component 
having  a  working  status  less  than  0.5  is  represented  by  a  hollow  symbol  to 
indicate  either  complete  or  partial  loss  of  functionality.  As  expected,  all 
electricity  transformers  are  shut  down,  the  electricity  network  is  disabled, 
while  only  local  electricity  generators  can  work  based  on  fuel,  providing 
limited  power  supply  to  nearby  communities.  Furthermore,  the  water  net¬ 
work  is  also  disrupted,  while  only  local  commercial  water  vendors  that  can 
pump  water  from  wells  using  generators  still  function.  As  a  result,  educa¬ 
tional  institutions  and  healthcare  facilities  are  all  disrupted  due  to  the 
shortage  of  power  and  water. 

The  social  impacts  are  estimated  as  follows.  In  total,  62.4%  of  the  water 
nodes,  0.0%  of  the  food  nodes,  0.0%  of  the  fuel  nodes,  and  100.0%  of  edu¬ 
cation  and  healthcare  nodes  are  disrupted  (i.e.,  A"ater  =  17.5 ,  Afod  =  0.0 , 
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Afel  =  0.0 ,  A“lucallon  =  84.0 ,  and  A'’calll'care  =  4.0 ).  There  is  no  demand  loss  for 
water  and  food,  but  all  demand  for  education  and  healthcare  are  lost.  In¬ 
terestingly,  after  disruption,  the  average  access  cost  for  food  even  de¬ 
creases  from  0.35  to  0.32.  This  is  mainly  due  to  the  loss  of  education  and 
healthcare  needs,  which  relieves  traffic  congestion  in  some  parts  of  the 
city.  The  overall  average  access  cost  for  water  increases  from  0.25  to  1.42, 
and  4.3%  of  population  lost  access  to  water,  which  is  caused  by  the  cost  in¬ 
crease  of  queuing  at  remaining  water  facilities.  The  spatial  distribution  of 
changes  in  communities’  water  access  costs  shows  that  costs  increased  sig¬ 
nificantly  around  several  of  the  commercial  water  providers  as  competi¬ 
tion  increases  for  those  limited  resources. 


Figure  14.  Water  access  cost  ratios  for  Beta  Model,  Scenario  0. 


Water  access  cost  ratio  after  vs.  before  disruption: 
Scenario  0 


Explanation 

■>T  Electricity-  substations 
e  Electricity-  generators 
e  Electricity- transmission  pylons 
Electricity-  transformers 
Fuel-  central  depot 
Fuel-  tanks 

Water-  treatment  plants 
Water-  towers 
Water-  tanks 
Water-  vendor 
Food-  processing  facilities 
Food-  markets 

- Major  roads 

—  Streams 
[  Water  bodies 

Water  Access  Cost  Ratio 

High  :  1411.63 

Low  :  1.09205 


after  cascading  disruption 


3.3.2  Beta  model  sensitivity  analysis 

The  above  case  study  (Bo)  is  denoted  as  benchmark  Scenario  o,  and  now 
six  additional,  representative  scenarios  are  examined  under  different  pa¬ 
rameter  settings,  focusing  on  three  societal  impacts  in  the  food  and  water 
sectors  and  using  the  Beta  Model.  Results  of  all  seven  scenarios  are  sum¬ 
marized  in  Table  3. 
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In  Scenario  l,  the  consideration  of  queuing  cost  at  resource  location  was 
removed.  Since  the  transportation  cost  is  the  total  access  cost,  the  reduc¬ 
tion  of  traffic  leads  to  a  more  noticeable  decrease  in  the  food  access  cost. 
There  is  no  significant  change  in  the  water  access  cost  after  disruption, 
which  means  the  negative  effect  of  water  facilities  failure  is  cancelled  out 
by  the  reduced  traffic  congestion. 


Table  3.  Societal  impacts  in  the  food  and  water  sectors  under  different  scenarios, 

using  the  Beta  Model. 


Scenarios 

Food 

Water 

Access 

cost 

increase 

Demand 

loss 

#  Nodes 

failed 
(total  11) 

Access 

cost 

increase 

Demand 

loss 

#Node 
failed  (total 
28) 

0:  Case  Study 

-7% 

0.0% 

0.0 

458% 

4.3% 

17.5 

1:  No  Queuing 

-36% 

0.0% 

0.0 

0.5% 

0.0% 

18.0 

2:  Kcap  =5-0 

-49% 

0.0% 

0.0 

-25% 

0.0% 

17.5 

3:  Kcap  =10.0 

-54% 

0.0% 

0.0 

-35% 

0.0% 

17.5 

4:  Initial  Water 

-27% 

0.0% 

0.0 

20% 

2.9% 

9.3 

5:  Initial  Fuel 

22% 

1.3% 

0.5 

157% 

5.9% 

10.6 

6:  Water  and  Fuel 

-2% 

0.0% 

0.3 

860% 

8.7% 

16.2 

Scenarios  2  and  3  examine  the  impact  of  facility  maximum  resource  ca¬ 
pacity  Xj .  It  is  expected  that  an  infrastructure  system  with  higher  maxi¬ 
mum  resource  capacity  would  have  greater  resiliency  to  system  disruption. 
In  Scenario  o,  let  the  resource  users  visit  their  nearest  resource  locations 
and  compute  the  resource  capacity  of  the  resource  locations  based  on  the 
demand  by  multiplying  a  factor  tccap  of  1.1.  In  Scenario  2  and  3,  Kcap  is 

changed  to  5.0  and  10.0,  respectively,  to  consider  an  infrastructure  system 
with  moderate  and  high  resource  capacity.  Similar  to  Scenario  1,  as  the 
queuing  cost  is  reduced,  the  transportation  cost  takes  a  greater  portion  in 
the  overall  access  cost;  hence,  the  alleviation  of  traffic  congestion  results  in 
a  greater  percentage  reduction  in  access  cost.  Decreases  of  49%  and  25% 
respectively  in  food  and  water  access  cost  were  observed  in  Scenario  2,  and 
54%  and  35%  decreases,  respectively,  were  observed  in  Scenario  3. 


Scenario  4  explores  a  different  initial  disruption  at  a  water  treatment 
plant.  Counterintuitively,  although  the  water  treatment  plant  is  the  key  in- 
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frastructure  in  the  water  layer,  shutting  it  down  does  not  significantly  af¬ 
fect  the  entire  system.  Although  4.0  water  towers  (due  to  support  failure), 
3.0  water  tanks,  and  2.3  water  vendors  (due  to  resource  failure)  are  dis¬ 
rupted,  the  water  access  cost  only  slightly  increases  (by  20%),  and  there  is 
a  2.9%  water  demand  loss.  These  numbers  are  mainly  because  of  the  high 
inelasticity  of  resource  demand  of  facilities  makes  resource  failure  very  un¬ 
likely  to  happen.  Further,  no  disruptions  in  the  food  layer  are  found,  and 
there  is  a  decrease  of  27%  in  food  procurement  cost  due  to  traffic  conges¬ 
tion  that  causes  increased  travel. 

In  Scenario  5,  the  initial  disruption  is  assumed  to  occur  at  the  central 
fuel  depot,  which  directly  cuts  down  the  fuel  supply  to  the  entire  system. 

As  a  result,  ten  local  water  vendors  are  completely  disrupted  (due  to  re¬ 
source  failure),  and  the  status  of  several  water  tanks  are  impaired  due  to 
increased  congestion.  Similarly,  the  food  markets  are  also  affected  due  to 
congestion,  with  a  total  status  reduction  of  0.5.  The  community  can  still 
survive  well  without  fuel  supply,  whereas  the  increases  in  access  cost  for 
water  and  food  have  a  moderate  effect.  Therefore,  the  fuel  layer  does  not 
play  a  critical  role  in  the  infrastructure  system  if  the  power  grid  and  water 
network  still  function  well. 

In  Scenario  6,  the  water  treatment  plant  and  the  fuel  depot  are  both 
failed  in  the  initial  disruption.  The  result  turns  out  to  be  more  devastating 
than  the  simple  summation  of  the  results  where  the  two  facilities  are  dis¬ 
rupted.  The  water  access  cost  is  increased  by  860%,  with  8.7%  of  the  popu¬ 
lation  having  lost  water  access.  The  demand  loss  appears  to  be  low,  while 
there  is  a  dramatic  increase  in  access  cost,  because  of  both  the  inelasticity 
of  resource  demand  and  a  relatively  insufficient  resource  supply  in  the  in¬ 
frastructure  system. 
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4  Conclusions 

The  goal  of  this  research  was  to  develop  a  system  that  can  be  used  by  mili¬ 
tary  planners  and  strategists  to  incorporate  longer-term  thinking  into  con¬ 
flict  analysis,  so  that  both  humanitarian  assistance  and  reconstruction 
strategies  can  be  better  visualized  during  earlier  phases  of  the  military 
planning  process. 

In  this  report,  a  holistic  mathematical  model  has  been  proposed  to  evalu¬ 
ate  the  vulnerability  of  an  urban  infrastructure  system  against  the  threats 
of  cascading  failures.  The  complex,  interdependent  infrastructure  systems 
are  modeled  as  a  multi-layered  network,  where  each  functioning  infra¬ 
structure  unit  is  modeled  as  a  node  in  a  layer  for  its  type.  Two  types  of  in¬ 
frastructure  failure  mechanisms  are  modeled  to  estimate  the  cascading 
failure  under  a  disaster:  support  failure  (e.g.,  disruption  due  to  loss  of  di¬ 
rect  physical  connections)  and  resource  failure  (e.g.,  disruption  due  to  loss 
of  resource  supply).  To  this  end,  a  revised  network  equilibrium  model  that 
incorporates  queuing  and  congestion  is  formulated,  and  the  ultimate  sys¬ 
tem  state  is  obtained  regarding  infrastructure  failure  and  community  be¬ 
havior  is  obtained. 

Through  a  case  study  done  for  this  work  and  based  on  information  from 
Maiduguri,  Nigeria,  several  interesting  insights  were  obtained.  First,  a  sys¬ 
tem  with  greater  resource  capacity  is  more  resilient  to  disruptions  because 
the  effect  of  resource  shortage  is  mitigated.  Second,  isolated  disruption 
happening  at  some  “seemingly”  critical  infrastructures  (such  as  an  up¬ 
stream  water  treatment  plant)  may  not  severely  affect  the  entire  system, 
because  the  city  may  be  equipped  with  a  considerable  amount  of  inde¬ 
pendent  backup  supplies  (e.g.,  commercially-run  water  wells).  Third,  un¬ 
der  some  scenarios,  maintaining  the  functionality  of  some  infrastructures 
that  are  struggling  to  survive  may  not  benefit  the  society.  Shutting  them 
down  may  help  relieve  resource  competition  and  guarantee  supply  of  nec¬ 
essary  resources  to  other  parts  of  the  system.  In  general,  mutual  interac¬ 
tions  among  infrastructures  and  communities  lead  to  high -vulnerability 
risks  to  the  system,  especially  when  the  initial  disruption  is  difficult  to  pre¬ 
dict  or  evaluate  in  advance.  This  finding  emphasizes  the  importance  of  us¬ 
ing  a  holistic  model  to  correctly  capture  the  interdependent  relationships 
between  the  physical  infrastructure  systems  and  the  population,  when 
evaluating  the  system  reliability  and  societal  impacts. 
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Future  research  could  be  conducted  in  a  number  of  directions.  First,  the 
model  can  serve  as  a  building  block  to  identify  the  most  critical  infrastruc¬ 
ture  in  a  given  system.  Based  on  this,  geospatially-informed  reinforcement 
or  interdiction  strategies  can  be  sought  that  will  best  protect  an  urban  sys¬ 
tem  as  a  whole  or  at  least  certain  vulnerable  communities  within  it.  Sec¬ 
ond,  the  model’s  framework  only  uses  modeled  values  of  time  to  measure  a 
community’s  generalized  cost;  this  can  be  naturally  expanded  by  consider¬ 
ing  additional  cost  components  (e.g.,  price  of  resources).  In  addition,  the 
model’s  framework  is  deterministic;  a  stochastic  infrastructure  failure 
scheme,  possibly  based  on  Markovian  transitions,  may  help  reveal  more 
insightful  and  realistic  results.  Finally,  for  effective  use  of  this  type  of 
model  for  planning  purposes,  it  will  be  critical  to  better  understand  the 
sensitivity  of  model  outputs  to  the  level  of  detail  and  structural  differences 
in  the  input  data  that  define  the  joint  infrastructure-community  network. 
Close  collaboration  between  potential  end-users  in  the  military  planning 
community  and  model  developers  is  recommended  as  model  development 
continues. 
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